IBBHiBWi 


■01'%';'; 


M 11111111 11111  i 

''.'•  .--.i. '■•■■'■■    ''■"■.'■  -■•■.■■■.•;■  '■'1  ■■■■;;■  ■ 

,■■.■■'-,..■  :■',-:  ''';!■::.■■■■'■'■'■■.'■■,■' ■■''•■^. 
•■'',""■■■,■'■'  •■■.■■A->'  vy. mm  .v.... i&Hjr  /.'■!■;:■  ■ 


CuWUU 


;*SM'.(lf> 


HilliEi 


■^r 


'm'K 


"(:■/;':! 


,r  ]./.:. 


mm 


V  '  ■"■  DKXwamQWO -■■;«•  ".■■■:■■* ■-■!.' :.■■'•,■,!;■ 
.■V:^;[^^:.'V.1iw-';>!'.^k,,:^:;;>\,!:,::'.''V  ,.' 


Dudley  Knox  Library,  NPS 
Monterey,  CA  93943 


NPS68-83-003 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

A 

STATISTICAL  APPROACH  FOR   DETERMINING 

SUBSURFACE  THERMAL  STRUCTURE 

FROM  SEA   SURFACE  TEMPERATURE 

IN  THE   NORTHEAST    PACIFIC    OCEAN 

by 

Terry  A.    Howell 

June    1983 

Th< 

3SiS 

Advisor:                                      G.    H.    Jung 

Approved  for  public  release;  distribution  unlimited 

Prepared  for:   Naval  Ocean  Systems  Center 

Code  531 

San  Diego,  California  92152 


T20S0 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  John  J.  Skelund,  USN         David  A.  Schrady 
Superintendent  Provost 


This  thesis  was  prepared  in  conjunction  with  researcn 
supported  in  part  by  Naval  Ocean  Systems  Center  under  work 
order  N6600183WR0001*3, 

Reproduction  of  all  or  part  of  this  report  is  not  authorized 
without  permission  of  the  Naval  Postgraduate  School. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.    REPORT  NUMBER         ~ 

NPS68-83-003 


2.  GOVT  ACCESSION  NO. 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.    TITLE  (and  Subtitle) 

A  Statistical  Approach  For  Determining 
Subsurface  Thermal  Structure  From  Sea 
Surface  Temperature  In  The  Northeast 
Pacific  Ocean 


5.     TYPE  OF   REPORT  &   PERIOD  COVERED 


Master's  Thesis 
June  1982 


PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR*; 

Terry  A.  Kowell  in  conjunction  with 
Glen  H.  Jung  and  Calvin  R.  Dunlap 


B.     CONTRACT  OR  GRANT  NUMBERf*.) 


i      PERFORMING  ORGANIZATION  NAME  AND  AOORE5S 

Naval  Postgraduate  School 
Monterey,  California  939^0 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

N6600t83WR0001+3 


tl.     CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Ocean  Systems  Center 
Code  531 

San  DiPgo,  California 2211 


12.     REPORT  DATE 

June    1983 


13.     NUMBER  OF  PAGES 


14.    MONITORING- AGENCY  NAME  *   AOORESSfi/  different  from  Controlling  Office) 

CNO  Attn:  OP095 
Navy  Department 
Washington,  D.C.   20350 


193 


IS.     SECURITY  CLASS,   (ol  this  report) 

Unclassified 


1S«.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.     DISTRIBUTION  STATEMENT  (ol  thla  Report) 

Approved  for  public  release;  distribution  unlimited 


t7.     DISTRIBUTION  STATEMENT  (ol  the  abatraet  entered  In  Block  20.  II  dlllerent  /root  Report) 


U.    SUPPLEMENTARY  NOTES 


19.    KEY  WOROS  (Continue  on  rererae  aide  It  neeeeaary  and  identity  by  block  number) 

Subsurface  Thermal  Structure        Remote  Sensing 


Sea  Surface  Temperature 
Northeast  Pacific  Ocean 


Satellite  Oceanography 
Fronts 


20.     ABSTRACT  (  Continue  on  reveree  eldm  II  neceaaary  and  Identity  by  block  number) 

Bathythermograph  data  acquired  from  the  research  vessel  USN3 
SILAS  BENT  along  a  meridional  track  in  the  Northeast  Pacific 
during  September  1977  were  statistically  analyzed  to  determine 
possible  associations  between  the  subsurface  tnermal  structure 
and  sea  surface  temperature.   Strongly  correlated  variables 
(thermocline  gradients,  mixed  layer  depth,  and  locations  of  the 
seasonal  and  main  thermoclines)  within  the  vertical  temperature 


do  ,; 


FORM 
AN  73 


1473  EDITION  OF   I  NOV  85  IS  OBSOLETE 

S/N  0)02-  LF- 014-660) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF   THIS  PAGE  (When  Data  Snterec 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whun  Dot.  Bnfnd) 


#20.  ABSTRACT  (Continued) 

profile  were  used  in  linear  regression  methods  to  form 
empirical  relationships.  The  generated  equations  then  are 
utilized  to  define  the  subsurface  thermal  structure  from 
only  an  input  of  sea  surface  temperature.   Comparison  tests 
with  temporally  and  spatially  removed  BT  data  were  conducted 
with  results  indicating  successful  application  within  a 
water  mass  domain  with  uniformly  changing  characteristics. 


S'N  0102-  LF-  014-6601 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(TWi«n  Omtm  Bnturtd) 


Approved  for  public  release;  distribution  unlimited 

A  Statistical  Approach  For  Determining 
Subsurface  Thermal  Structure  From  Sea  Surface  Temperature 
in  the  Northeast  Pacific  Ocean 

by 

Terry  A.  Howell 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  197*+ 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  OCEANOGRAPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1983 


fi  2297 


Dudley  Knox  Library.  NPS 
Monterey.  CA  93943 


ABSTRACT 


Eathyt  heroic  graph  data  acquired  from  the  research  vessel 
USNS  SILAS  BENT  along  a  meridional  track  in  the  Northeast 
Pacific  during  September  1977  were  statistically  analyzed  to 
determine  possible  associations  between  the  subsurface  ther- 
mal structure  and  sea  surface  temperature.  Strongly  corre- 
lated variables  (ther mocline  gradients,  mixed  layer  depth, 
and  locaticns  cf  the  seasonal  and  main  thermoclines)  within 
the  vertical  temperature  profile  were  used  in  linear  regres- 
sion methcds  tc  fcrm  empirical  relationships.  The  generated 
equations  then  are  utilized  to  define  the  subsurface  thermal 
structure  from  only  an  input  of  sea  surface  temperature. 
Comparison  -sets  with  tenpcrally  and  spatially  removed  BT 
data  were  ccnducted  with  results  indicating  successful 
application  within  a  water  mass  domain  with  uniformly  chang- 
ing   characteristics. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 15 

A.  IMPORTANCE  OF  DETERMINING  OCEAN  THERMAL 
STRUCTURE 15 

B.  SATELLITE  DETER KINATIO N  OF  SUBSURFACE  THERMAL 
STRUCTURE 16 

C.  PURPCSE  AND  GOAIS  OF  THESIS 17 

II.  GCEANCGBAEHIC  CHARACTERISTICS  OF  THE  NORTHEAST 

EACIFIC  OCEAN   19 

A.  WINE  SYSTEMS  ANC  CURRENTS 19 

B.  WATER  MASSES  ANE  FRONTS 21 

1.  SUEARCTIC  WATER  MASS  CHARACTERISTICS   ...  21 

2.  SUETROPIC  WATER  MASS  CHARACTERISTICS   ...  25 

3.  TRANSITION  ZCNE 28 

4.  DENSITY  STRUCTURE  AND  THE  SUBARCTIC  FRONT  .  29 

III.  THERMAL  STRUCTURE 32 

A.  THESKAL  STRUCTURE  IN  THE  NORTHEAST  PACIFIC 

OCEAN 32 

B.  INTERNAL  WAVE  EJEECTS  ON  THERMAL  STRUCTURE   .  .  36 

C.  HORIZONTAL  SURFACE  THERMAL  STRUCTURE   38 

IV.  EXPERIMENT 40 

A.  OBJECTIVE 40 

B.  LATA  ACQUISITION 42 


C.  ASSESSMENT    OF    TEMPERATURE    STRUCTURE    45 

D.  STATISTICAL    PROCEDURE    61 

E.  REGRESSION    EQUATION    ANALYSIS       81 

F.  THERMAL   PROFILE    TEST    AND    COMPARISON    ANALYSIS       .     87 

1.  TEST    1 87 

2.  TEST    2 90 

2.       TEST    3 10  1 

G.  REVIEW    CF    PREDICTED    THERMAL    PROFILE    RESULTS     .        107 
V.             CONCLUSIONS   AND    RECCMMENDA TIONS    113 

A.  CONCLUSIONS 113 

B.  RECOMMENDATIONS 115 

APPENDIX    A.       THERMAL    VARIAELES    VERSUS    LATITUDE/DISTANCE 

ALONG    TRACK 118 

APPENDIX    E.       EICMEDICAL    COMPUTER    PROGRAMS     (DIXON    AND 

EBCWN,     1979)     132 

APPENDIX    C.       REGIONAL    SCATTER    DIAGRAMS    -    S3T    VERSUS 

THERMAL    VARIAELES       134 

APPENDIX    D.       REGIONAL    SCATIER    DIAGRAMS    -    BOT    VERSUS 

THERMAL    VARIAELES       147 

APPENDIX    E.       REGIONAL    SCATTER    DIAGRAMS    -    BUZ    VERSUS 

THEFMAL    VARIAELES       158 


APPENDIX     F.       OBSERVED    AND    PREDICTED    THERMAL    PROFILES 


AIONG   TRACK 


168 


LIST    OF    REFERENCES 


188 


INITIAL  DISTRIEOTION  LIST 


191 


LIST   OF    TABLES 


TABLE   I.  locations   and   Times    of    XBT   Measurements      .    .    51 

TABLE  II.  Values    of   SSI,    MLDr    Thermoclina   Gradients 

and    T460 52 

TABLE   III.         Values    of  Th€rmocline   Temperatures   and 

Depths 53 

TABLE   IV.  Correlation    Coefficients   of   Thermal 

Variables 64 

TABLE   V.  Best-Pit   Correlation   Coefficiants   of   Thermal 

Variables 79 


Figure  1. 

Figure  2. 

Figure  3. 

Figure  4. 

Figure  5. 

Figure  6. 

Figure  7. 

Figure  8. 

Figure  9. 

Figure  10. 

Figure  1 1 . 

Figure  12. 

Figure  13. 

Figure  14. 

Figure  15. 

Figure  16. 

Figure  17. 

Figure  18. 


LI  SI    OF    FIGURES 

Location    of    OfcS-P   and    North  Pacific   Ocean 
Surface    Currerts   20 

Characteristic  Temperature   and   Salinity 
Profiles 22 

Upper  Zone  Domains 26 

Lower  Zone   Domains 27 

Growth   and   Decay   of    the   Thermccline    at   Ocean 
Station    "P"       34 

Typical    Temoerature    Structures    at    Ocean 
Station    "P* 37 

Surface    Mean    Temoerature   Structure    For 
September    1977    ". 39 

Lccation    of    XET   DATA    During   8-11    September 
1977 44 

Meridional  Cress-section   of  Temperature  at 
145W 46 

Model  of    Vertical   Temperature    Structure.       .    .    49 

Pacific    Ocean    Sea    Surface   Temperature   on    10 
September    1977 55 

Ccaparison   of    ET5,    DT15,    DT50    and   DT100 
Versus    Latitude 58 

Comparison   of    MLD,    BOT    and   BUZ    Versus 
Latitude 6  0 

Ccmoarison   of    SST r    TBOT    and  TBUZ    Versus 
Latitude 66 

Ccipariscn   of    $LD   and    SST    Versus    Latitude.       .  67 

Comparison   of    EOT    and    SST    Versus   Latitude.       .  69 

Ccirparison   of   EUZ   and   SST    Versus   Latitude.       .  69 

Division    cf    Regions    Along    Track 71 

9 


Figure  19.   Regional  Scatter  Diagram  -  DZ  Versus  30T.  .  .  77 

Figure  20.   Regional  Scatter  Diagram  -  DZBUZ  Versus  BUZ.   78 

Figure  21.   Observed  and  Predicted  Thermal  Profiles  -  X3T 

Station  339 88 

Figure  22.   Location  of  Original  (•)  and  Test  2  (X)  XBT 

Drccs 92 

Figure   23.      Observed    and    Predicted   Thermal    Profiles  at 

41.  5 N,     140.3 W 94 

Figure   24.      Observed    and    Predicted    Thermal    Profiles   at 

42. 3N,     150. 0W 95 

Figure   25.      Observed    and    Predicted   Thermal    Profiles   at 

41. 9N,     148. 7W      97 

Figure   26.      Observed    and    Predicted   Thermal   Profiles   at 

42. 3N,     144. 0W 93 

Figure   27.      Observed    and    Predicted   Thermal   Profiles  at 

41.  1N,     143. 2H „ 99 

Figure   28.      Observed    and    Predicted   Thermal   Profiles   at 

CKS-P 100 

Figure   29.      Lcca-ion    of    Original    (•)    and   Test    3     (X)    XBT 

Data 102 

Figure  30.   Observed  and  Predicted  Thermal  Profiles  for 

182030Z  SEP  77 104 

Figure  31.   Observed  and  Predicted  Thermal  Profiles  for 

240000Z  SEP  77 106 

Figure  32.   Observed  and  Predicted  Thermal  Profiles  for 

300000Z  SEP  77 109 

Figure  33.   Example  Predicted  Profile  Error  Due  to 

"afternoon  Effect" 111 

Figure   A.1.    SST  Versus  Latitude/Distance 118 

Figure   A.  2.    HLD  Versus  Latitude/Distance 119 

Figure   A. 3.    DT5  Versus  Latitude/Distance 120 


10 


Figure      A. 4.         DT15    Versus   La-it ude/Distance 121 

Figure      A. 5.        DT50    Versus   Latitude/Distance 122 

Figure      A. 6-        DT100    Versus  Latitude/Distance 123 

Figure      A.  7.         BCT  Versus    Latitude/Distance 124 

Figure      A. 8.        TEOT    Versus  Latitude/Distance 125 

Figure      A.  5.        DT  Versus    Latitude/Distance 126 

Figure      A.  10.         EZ   Versus    Latitude/Distance 127 

Figure      A.  11.         BUZ    Versus   Latitude/Distance. 128 

Figure      A.  12.        TBUZ    Versus   Latitude/Distance 129 

Figure      A.  13.         ETBOZ   Versus   Lat  itude/Dist  ancs 130 

Figure      A.  14.         DZBCJZ   Versus   Latitude/Distance 131 

Figure     C.1.        Regional    Scatter    Diagram   -    SST   Versus 

ELD 134 

Figure  C.2.    Regional  Scatter  Diagram  -  SST  Vsrsus 

CT5 135 

Figure   C.3.    Reqional  Scatter  Diaaram  -  SST  Versus 

D115. " 136 

Figure     C.4.         Regional    Scatter    Diagram  -    SST   Versus 

DT5C 137 

Figure      C.5.         Reqional    Scatter    Diagram  -    SST   Versus 

DT100 138 

Figure  C.6.    Regional  Scatter  Diagram  -  SST  Versus* 

ECT 139 

Figure  C.7.    Regional  Scatter  Diagram  -  SST  Versus 

1EOT 140 

Figure      C.8.         Regional    Scatter    Diagram  -    SST    Versus    DT.  141 

Figure      C.9.         Regional    Scatter    Diagram  -    SSI    Versus    DZ.  142 


11 


Fiaure      C.10.         Regional    Scatter    Diagram   -    SST    Versus 

EUZ. 143 

Figure      C.11.         Eegional    Scatter    Diagram   -    SST    Versus 

TEDZ 144 

Figure     C.12.        Regional    Scatter    Diagram    -   SST    Versus 

CTEUZ 145 

Figure     C.13.         Regional    Scatter    Diagram   -    SST    Versus 

DZBUZ 146 

Fiaure      D.1.         Regional    Scatter    Diagram  -    30T   Versus 

DT5 147 

Figure      D.2.         Regional    Scatter    Diagram  -    30T    Versus 

E115 148 

Figure   D.3.    Regional  Scatter  Diagram  -  BOT  Versus 

CT5C 149 

Figure      D.4.        Regional   Scatter    Diagram   -    BOT    Versus 

ET100 150 

Figure   D.5.    Regional  Scatter  Diagram  -  BOT  Versus 

TEOT. 151 

Figure   D.6.    Regional  Scatter  Diagram  -  BOT  Versus  DT.  152 

Figure   D.7.    Regional  Scatter  Diagram  -  BOT  Versus  DZ.  153 

Figure   D.8.    Regional  Scatter  Diagram  -  BOT  Versus 

EUZ 154 

Figure      D.9.         Regional    Scatter    Diagram  -    30T   Versus 

TEUZ 155 

Figure      D.10.         Eegional    Scatter    Diagram   -    BOT    Versus 

D130Z 156 

Figure      D.11.         Eegional    Scatter    Diagram   -    BOT    Versus 

EZBUZ 157 

Figure   E.1.    Regional  Scatter  Diagram  -  BUZ  Versus 

ET5 153 

Figure   E.2.    Regional  Scatter  Diagram  -  BUZ  Versus 

DT15 159 


12 


Figure      E.3.        Regional    Scatter    Diagram  -    BUZ   Versus 

ET50. 160 

Figure   E.4.    Regional  Scatter  Diagram  -  BUZ  Versus 

ET100 161 

Figure      E.5.         Reaional    Scatter    Diagram  -    BUZ    Versus 

ISO!: 162 

Figure      E.6.        Regional    Scatter    Diagram  -    BUZ    Versus    DT.     163 

Figure      E.7.         Regional    Scatter    Diagram  -    BUZ    Versus    DZ.    164 

Figure      E.8.        Regional    Scatter    Diagram  -    BUZ   Versus 

TEUZ 165 

Figure      E.9.         Regional   Scatter    Diagram  -    BUZ    Versus 

DTBUZ 166 

Figure      E.  10.         Regional    Scatter    Diagram   -    BUZ    Versus 

DZBUZ 167 

Figure      F.1.         Observed    and  Predicted   Thermal    Profiles    at 

XET    Sta.    340 168 

Figure      F.2.        Observed    and   Predicted   Thermal   Profiles    at 

XET  Sta.  341 169 

Figure   F.3.    Observed  and  Predicted  Thermal  Profiles  at 

XET    Sta.    342 170 

Figure      F.4.        Observed    and  Predicted   Thermal   Profiles   at 

XET  Sta.  343 171 

Figure   F.5.    Observed  and  Predicted  Thermal  Profiles  at 

XET  Sta.  344 172 

Figure  F.6.    Observed  and  Predicted  Thermal  Profiles  at 

XET  Sta.  345 173 

Figure   F.7.    Observed  ar.d  Predicted  Thermal  Profiles  at 

XET  Sta.  34  6 174 

Figure   F.8.    Observed  and  Predicted  Thermal  Profiles  at 

XET  Sta.  347 175 

Figure   F.9.    Observed  and  Predicted  Thermal  Profiles  at 

XET  Sta.  348 176 


13 


Figure      F.10.         Observed    and   Predicted  Thermal   Profiles 

at    XBT    Sta.     349 177 

Figure      F.11.        Observed    and   Predicted  Thermal   Profiles 

at  XET  Sta.  350 173 

Figure  F.12.    Observed  and  Predicted  Thermal  Profiles 

ax  XET  Sta.  351 179 

Figure  F.13.    Observed  and  Predicted  Thermal  Profiles 

at    XBT    Sta.     352 180 

Figure      F.14.        Observed    and   Predicted  Thermal   Profiles 

at    XBT    Sta.     353 131 

Figure      F.15.        Observed    and   Predicted  Thermal   Profiles 

at  XET  Sta.  354 182 

Figure  F.16.    Observed  and  Predicted  Thermal  Profiles 

at  XET  Sta.  355 183 

Figure  F.17.    Observed  and  Predicted  Thermal  Profiles 

at    XBT    Sta.     356 134 

Figure      F.18.         Observed    and   Predicted  Thermal   Profiles 

at    XBT    S^a.     357 185 

Figure      F.19.        Observed    and   Predicted   Thermal    Profiles 

at    XET    Sta.     358 186 

Figure      F.20.         Observed    and   Predicted  .Thermal    Profiles 

at    XET    S^-.a.     359 187 


14 


I.        INTRODUCTION 

A.        IMPORTANCE    OF   DETERMINING    OCEAN    THERMAL    STRUCTURE 

An  important  physical  property  and  product  cf  various 
dynamical  processes  occurring  in  -he  oceans  is  the  thermal 
structure.  Although  it  is  derived  from  only  a  few  mecha- 
nisms such  as  air-sea  thermal  energy  exchange,  forced  and 
free  convective  mixing  and  horizontal  ocean  thermal  advec- 
tion,  the  interaction  of  these  processes  along  with  ether 
hydrcdynamic  motions  can  become  immensely  complicated. 
Since  the  temperature  profile  has  important  implications  for 
climatic  weather  prediction,  commercial  fisheries  and  mili- 
tary acoustic  surveillance,  to  men-ion  a  few,  the  ability 
to    predict   this    structure    is   of    major   importance. 

Climatological  atlases  present  an  averaged  profile  for  a 
certain  period  and  location.  However,  this  results  ir.  a 
generalized  tieriral  structure  which  is  not  accurate  enough 
to  accomplish  the  precise  needs  cf  underwater  sound  propaga- 
tion forecasting.  Recent  developments  in  numerical  modeling 
of  the  upper-ocean  temperatures  have  produced  vast  improve- 
ments in  the  analysis  and  prediction  of  mixed  layers  and 
other  synoptic  oceanic  features  (Clancy  et  al.  ,  1981).  It 
is    inevitable   that   this    method   will   reach    the    sophistication 


15 


cf  today's  meteorological  forecasting,  but  with  the  complex- 
ity cf  the  physical  processes  involved  and  with  present  data 
sparsity  ir.  many  ocean  areas,  many  inaccuracies  still  exist. 
The  bathythermograph  (BT)  provides  an  accurate 
temperature  profile  at  a  point  location  in  the  upper  ocean. 
Since  the  vertical  thermal  structure  over  a  large  area  pro- 
vides ths  main  factor  which  affects  the  propagation  charac- 
teristics cf  underwater  acoustics,  the  need  for  many 
temperature  profiles  in  Navy  ASW  operations  is  obvious.  The 
expense  involved  in  ship  and  aircraft  deployment  of  BT's  for 
an  ocean  area  of  interest  often  prohibits  the  use  of  this 
most  accurate  means  of  obtaining  a  detailed  thermal 
description. 

B.        SATELLITE    DETERMINATION    OF    SUBSURFACE    THERMAL    STRUCTURE 

The  study  cf  the  oceans  by  satellites  nas  become  a  major 
arena  for  scientific  scrutiny  and  investigation.  Already  it 
has  given  added  insight  into  the  numerous  processes  that 
govern  oceanic  behavior.  The  Gulf  Stream  eddies  and  fronts, 
the  coastal  upweiling  off  the  Somali  and  Arabian  coasts,  and 
other  applications  in  biological,  geophysical  and  dynamical 
fields  are  becoming  well  known.  However,  in  spite  of  great 
progress,  satellite  remote  sensing  techniques  for  observing 
oceanic   processes  are   far    from   being   fully   realized. 
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The  fundamental  problem  which  continuously  arises, 
besides  atmospheric  effects  which  will  not  be  discussed,  is 
that  satellite  observations  are  made  at  the  surface  alcne 
and  that  all  conclusions  about  subsurface  processes  are 
obtained  from  inference  rather  than  by  direct  measurement. 
Consequently,  the  applications  of  satellites  for  oceanogra- 
phy are  limited  to  those  phenomena  that  generate  surface 
signatures  detectable  by  electromagnetic  waves.  Once  these 
surface  manifestations  are  observed,  it  is  necessary  to  have 
knowledge  of  the  physical  processes  involved  to  interpret 
their  meaning.  Currently  there   are    few      methods   available 

where  a  qualitative  "feel"  for  subsurface  features  can  be 
directly  expressed  as  empirical  relationships  from  remotely 
sensed  measurements  at  the  sea  surface.  If  a  linkage  is 
found  to  exist  between  the  air-sea  interface  and  the  thermal 
structure  below,  a  satellite  would  then  be  the  perfect  plat- 
form to  collect  sea  surface  data  rapidly  over  all  the 
world's  oceans.  This  surface  information  could  then  be  used 
to    infer   the    subsurface    vertical    temperature    structure. 

C.        EUHPOSE    AKE    GCALS    OF    THESIS 

The  primary  purpose  fcr  the  ongoing  research  into  the 
prediction   of   subsurface    thermal      structure   is   to    understand 
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the  processes  and  relationships  that  govern  the  ocean's 
behavior.  Bcles  that  are  performed  by  ships  and  aircraft 
(via      expendable     bathythermographs)  in      determining      the 

temperature  profiles  could  be  played  better  by  satellites, 
provided  a  "link"  is  found  to  read  the  below-surf ace  events 
from  the  sea  surface  signatures.  Therefore,  if  strong  cor- 
relations and  relationships  could  be  ascartained  between 
satellite-derived  sea  surface  temperatures  and  vertical 
temperature  profiles,  then  the  areas  of  acoustical  oceanog- 
raphy and  naval  tactical  applications  would  benefit  greatly. 
The  goal  of  this  thesis  is  to  determine  the  possible 
correlations  that  exist  between  the  subsurface  thermal 
structure  and  the  temperature  at  the  surface.  Then,  from 
these  relationships  and  use  of  regression  formulas, 
predicted  vertical  temperature  profiles  could  be  computed 
from    an    input    of    sea    surface   temperature   alone. 
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II.   CCEANOGEAPHIC  CHARACTERISTICS  OF  THE 


A.   RIND  SYSTEMS  AND  CURRENTS 

The  Northeast  Pacific  has  been  studied  extensively  for 
many  years  and  it  has  revealed  a  complex  arrangement  of 
environmental  conditions.  Severe  storms,  high  winds,  varied 
water  masses,  numerous  currents  and  unique  temperature  and 
salinity  structures  are  a  normal  occurrence  for  the  region. 

In  the  winter  the  North  Pacific  is  usually  under  the 
influence  of  the  Aleutian  Low  in  the  atmosphere  while  for 
the  summer  months  a  high  pressure  system  dominates  (Uda, 
1963).  The  distribution  of  the  mean  surface  winds  are  con- 
trolled by  these  two  systems.  At  the  location  of  the  new 
retired  Ocean  Weather  Station  "P"  (OWS-P,  positioned  at  SON, 
145W  in  the  Gulf  of  Alaska)  the  prevailing  winds  are  from 
the  southwest  (Tafcata,  1565) .  From  the  influence  of  these 
wind  patterns,  the  surface  waters  flow  eastward  across  the 
Pacific  Ocean  in  the  vicinity  of  Station  "P".  This  flow 
consists  of  the  west  Wind  Drift  and  Subarctic  Current  which 
then  diverge  to  form  the  southward  flowing  California  Cur- 
rent  and  the  northward   flowing   Alaskan  Current   (Tabata, 
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1978) .    Fig.   1  shows  the  location  of  Ocean  Weather  Station 
"P"  and  the  major  surface  currents  of  the  North  Pacific. 


JSf< 


NOWTH  PACIFIC  OCEAN 

GENERAL  CIRCULATION 


Figure    1.      Location    of   OWS-P   and    North  Pacific   Ocean    Surface 
Currents    (From    Tully,    1960) . 
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B.        WATER    MASSES    AND    FRONTS 

The   Northeast      Pacific  can  be      divided    into      three    major 
regions   based      upon    their    water    characteristics.  Near   the 

vicinity  cf  CWS-P  and  northward  is  ths  Subarctic  Pacific 
Water  Mass  in  which  the  upper  100-150  m  layer  is  composed  of 
cold  and  lew-salinity  water.  To  the  soutn  warm  and  highly 
saline  water  exists  in  the  upper  layers  of  the  Subtropic 
Water  Mass.  In  fcetwean  and  along  the  coast  lies  a  region  in 
which  characteristics  from  both  water  masses  are  combined. 
This  is  known  as  the  Transition  Zone  where  there  is  a  marked 
gradation  of  cceanographic  properties  from  one  region  to 
another,  with  the  boundaries  separating  these  water  masses 
having  complicated  therrachaiine  structures  (Roden,  1970). 
The  location  cf  the  Subarctic  Front  between  the  two  water 
masses  remains  fairly  ccrstart,  from  UON  to  U5Nr  and  is 
largely   dependent  on    the    prevailing   winds     (Roden,     1975). 

1.       SUEARC1IC    WATER    WSS    CHARACTERISTICS 

The  irain  feature  cf  the  vertical  structure  cf 
oceanic  properties  is  the  existence  of  an  upper  zone,  a 
halccline  and  a  lower  zone  (Uda,  1963).  The  temperature, 
salini-y  and  density  profiles  of  these  three  zones  for  sum- 
mer   and   wirt=r   are  depicted   in   Fig.    2. 
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Figure   2.      Characteristic   Temperature    and    Salinity   Profile; 
(Frcm    Tully,     1964)  . 
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The  upper  zone,  which  extends  frcm  the  surface  tc  a 
depth  of  about  1C0  m,  generally  has  a  low  salinity  level  of 
less  than  33  g/kg  with  the  lowest  values  being  closest  tc 
the  coastal  recicns  in  the  late  spring  (Uda,  1963) .  In  the 
winter,  the  water  in  this  zone  is  thoroughly  mixed  with  all 
physical    properties    becoming   homogeneous.  During   the   sum- 

mer, due  tc  lighter  surface  winds,  a  shallow  isohaline  and 
isothermal  layer  exists  from  10  to  30  m  balow  the  sea  sur- 
face. Below  this  mixed  layer  is  a  strong  negative  thermo- 
cline  with  tenperature  decreasing  as  much  as  8C  in  20  m;  in 
the  same  zone  a  small  halccline  extends  above  another  near- 
ischaline  layer  to  the  fccttom  of  the  upper  zone  (Tabata, 
1965)  . 

Below  the  upper  zone  lies  the  main  haioclina  which 
is  permanent,  as  the  salinity  increases  by  as  much  as  1  g/kg 
from  100  tc  200  m  (Tabata,  1960).  The  temparature  generally 
decreases  with  depth  in  this  region,  but  at  rimes  a 
temperature  inversion  occurs  tc  a  depth  of  about  150  m. 
This  positive  temperature  gradient  is  in  a  stable  layer 
(density  increases  with  depth)  since  it  is  located  within 
the    main    halocline    (Uda,     1963). 
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In  the  lcwer  zone,  fcelow  the  halocline,  the  proper- 
ties of  the  water  change  gradually  with  depth.  Salinity 
increases  slightly  by  about  0.6  g/kg  in  1000  m  while  -the 
temperature  decreases  to  a  value  of  approximately  2 . 8C  at  a 
depth   of    1000    a    {Tabata,    1S65)  . 

The  variations  of  salinity  are  most  pronounced  in 
the  upper  zone  and  halocline  and  are  less  in  the  lower  zone. 
The  annual  salinity  variations  in  the  upper  zone  are  charac- 
terized by  a  decrease  in  the  upper  30  m  in  the  late  summer 
or  early  fall  with  a  iraxiirum  in  late  winter  or  early  spring. 
These  seasonal  fluctuations  result  primarily  from  the  evapo- 
ration and  precipitation  imbalances  common  to  the  region. 
Annual  variation  is  similar  at  great  depths,  but  lags  the 
surface   layers   in  time    (Tabata,     1965)  . 

At  the  surface,  several  features  are  observed  in  the 
salinity  pattern.  Basically,  the  trend  of  the  surface 
isohalines  appears  to  be  east-west  throughout  the  year  with 
major  changes  occurring  along  the  coast.  Also,  a  permanent 
salinity  maximum  of  about  33  g/kg  exists  near  the  center  cf 
the  Gulf  cf  Alaska,  south  of  Kcdiak  Island.  The  salinity 
decreases  southward  to  a  minimum  of  about  32.6  g/kg  in  the 
vicinity  of  EON  and  then  increases  further  south  into  the 
Subtropical   Hater   Mass    (Uda,    196  3). 


24 


The  "three- zone"  structure  of  the  water  column  is 
characteristic  of  the  deep  ocean  region  of  -he  Subarctic 
Domain.  However,  it  does  not  always  exist  in  near-coastal 
waters  as  current  systems,  coastal  winds  and  fresh  water 
runoff  frcm  land  can  modify  the  profiles.  As  shown  in  Figs. 
3  and  4,  this  coastal  regime  primarily  alters  the  upper 
zone.  Along  the  western  United  States  this  domain  exists 
out  tc  about  130W.  The  lcwer  zone  also  becomes  affected  by 
the  proximity  to  land  of  the  California  Undercurrent 
(Dcdimead    et    al.  ,    1962). 

2.       SUETRCEIC   WATER    MJSS    CHARACTERISTICS 

To  the  south  of  the  Subarctic  region  is  the  Sub- 
tropic  Mater  Mass  which  is  basically  warmer  and  more  saline 
than  the  Subarctic  Water.  In  the  upper  layer  the  high  val- 
ues cf  salinity  (greater  than  34  g/kg)  due  to  higher  evapo- 
ration decrease  to  a  minimum  at  200  to  600  m  and  are  typical 
of  the  values  throughout  the  year.  Sea  surface  temperatures 
remain  fairly  ccr.stant  with  an  isothermal  upper  zone  averag- 
ing about  20  m  in  thickness  below  the  sea  surface,  due  tc 
the  prevailing  winds.  Eelow  this  zone  exists  a  negative 
permanent  thermocline  reaching  dcwn  to  a  depth  of  500-600  m. 
Also,        ther=      is      no        characteristic     halociine     or      major 
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Figure   3.       Upper    Zone    Domains    (From    Dodimeal    et    al.,     1962). 
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Figure   4.       Lcwez    Zone   Domains    (From    Dodimead    et   al..    1962) 
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temperature  inversion  as  there  was  in  the  water  mass  to  the 
north  (Tabata,  1965).  The  deep  water  structure  is  basically 
continuous  and  similar  from  the  Subarctic  to  Subtrcpic 
regions    (Oda,    1963)  . 

3.       TRANSITION    ZONE 

The  Transition  Zens,  located  in  the  central  and 
eastern  North  Eacific  between  45N  and  35N ,  separates  Subarc- 
tic water  frcm  Subtropic  water.  It  consists  of  the  mixed 
water  originating  in  the  warm,  saline  Kurcshio  and  the  cold, 
low  salinity  Cyashio  currents  as  they  move  eastward  in  the 
West    Wind    Drift    (Uda,    1963). 

The  width  of  this  zonal  belt,  which  is  about  2  tc  4 
degrees  of  latitude,  and  the  strength  of  its  boundaries,  are 
primarily  determined  by  the  wind  stress  at  the  sea  surface. 
The  northern  boundary  is  characterized  by  many  temperature 
inversions  and  by  the  gradual  disappearance  of  the  Subarctic 
halocline.  At  the  southern  boundary  the  vertical  structure 
becomes  basically  isohaline  with  a  strong  negative  thermo- 
cline   in    the    upper   layer    (Seder.,     1970). 

The  Transition  Zone  is  primarily  a  property  of  the 
upper  ocean,  as  its  characteristics  fade  below  a  few  hundred 
meters.  At   the      surface     this      region   has      a      distinctive 
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meridional  (N-S)  temperature  gradient  of  approximately  1.0C 
per  100km  as  it  stretches  eastward  across  the  Pacific  Ocean. 
The  horizontal  surface  salinity  pattern  exhibits  similar 
structure  along  a  latitude  band  throughout  the  year  with 
values  increasing  as  much  as  0.8  g/kg  meridionally  (Uda, 
1963)  . 

The  upper  layers  are  often  characterized  by 
temperature  inversions  which  can  occur  because  of  the  exis- 
tence of  the  associated  halocline.  These  inversions 
decrease  in  magnitude  and  slowly  increase  in  depth  toward 
lower   latitudes. 

The  vertical  distribution  of  properties  at  interme- 
diate depths  shows  that  the  salinity  becomes  effectively 
isohaline,  with  a  slight  nirimum  developing  below  300  m  at 
the    southern    boundary    (Uda,   1963). 

u-       DENSITY    STRUCTURE    AND    THE    SUBARCTIC    FRONT 

The  density  cf  sea  water  is  calculated  primarily 
from  its  temperature  and  salinity  at  depth.  In  the  Subarc- 
tic waters,  salinity  governs  the  density  structure.  The 
seasonal  variability  illustrates  some  of  the  features  that 
occur  in  the  vertical  profiles.  In  winter,  basically  iso- 
pycnal  conditions      are   found     throughout   the      upper  zone      as 
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temperature  and  salinity  are  well  mixed.  In  summer,  a  minor 
pycnccline  corresponding  to  the  shallow  thermocline  is 
observed.  Hcwever,  a  large  and  permanent  pycnocline  occurs 
along  with  the  main  halocline.  In  the  lower  zone,  the  den- 
sity increases  gradually  with  depth  (Uda,  1963). 

In  Suttropic  Water,  the  density  is  controlled  by  the 
temperature.  Here  the  warm  surface  layer  of  about  10  to 
30  m  is  iscpycnal  frcm  the  mixing  of  the  prevailing  winds. 
A  permanent  pycnccline  exists  at  lower  depths  due  to  the 
thermocline  in  the  region  (Odar  1963)  . 

The  Transition  Region  has  a  complex  density  struc- 
ture as  numerous  minor  density  fronts  occur  at  depth  due  to 
the  temperature  and  salinity  inversions. 

The  Subarctic  Front  in  the  eastern  North  Pacific  is 
a  narrow,  meandering  band  located  in  the  latitude  range  of 
40N  to  U5N.  This  front  has  a  unique  structure  in  its  lack 
of  a  density  front  in  the  upper  100  m.  The  strong  horizon- 
tal temperature  and  salirity  gradients  are  in  almost  com- 
plete balance.  Below  this  upper  layer  a  density  front  of 
moderate  intensity  slopes  slightly  to  the  south  (Ecden, 
1975)  . 
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The  mixed  layer  depth  can  strongly  change  across  the 
Subarctic  Front.  Tc  the  north,  the  mixed  layer  extends  to 
the  top  cf  the  halocline  at  about  100  m.  To  the  south, 
there  is  nc  halocline  and,  during  the  stormy,  winter  season, 
the  mixed  layer  depth  can  extend  to  almost  300  m.  In  the 
summer  months  the  seasonal  thermocline  controls  the  depth  of 
the  mixed  layer  and  there  is  little  difference  across  the 
Subarctic  Front  (Eoden,  1S75). 
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III.       THERMAL    STRUCTURE 

A.        THERMAL    STRUCTURE    IN    THE    NORTHEAST    PACIFIC    OCEAN 

The  interaction  between  the  sea  and  atmosphere  is  the 
dominant  process  determining  temperature  structure  in  the 
upper  layers.  Below,  temperature  is  a  fairly  cons-ant  prop- 
erty with  the  structure  being  determined  by  internal  pro- 
cesses,   such   as   advection   and   mixing. 

Two  heating  cycles  occur  in  the  ocean:  diurnal  and  sea- 
sonal. In  the  daily  cycle  the  surface  layer  is  heated  dur- 
ing the  daylight  hours  and  cooled  during  the  night.  In  the 
higher  latitudes  during  the  heating  season,  i.e.,  the  middle 
six  months  of  the  year,  a  daily  net  heat  gain  occurs  (Uda, 
1963).  Fig.  5  illustrates  the  annual  cycle  of  heating  and 
cooling      and    its     relation   to      the      thermal   structure.  In 

March,  at  the  end  of  the  coding  season,  the  waters  are 
nearly  isothermal  to  the  top  of  the  haiociine  located  at 
about  100  m  in  depth.  With  the  beginning  of  the  heating 
seascn      in      April,  the      heat      is      accumulated        as      small 

temperature  inversions  are  gradually  mixed  downward  by  light 
winds.         During    subsequent   short    periods   of    high    winds    these 
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transients   are   quickly   eroded   and      a   shallow    isothermal    sur- 
face   layer   results. 

Eelow  the  surface  layer,  a  sharp  negative  thermocline  is 
formed  between  30  and  60  m  in  a  temperature  transition  layer 
adjacent    to      the    deeper    cold   waters.  From   the    2nd      cf    the 

heating   season    until    the    waters      are   again   nearly   isothermal 
in    March,    the   thermocline    deepens  to   about    100    m. 

In  this  ocean  region,  below  the  thermocline  there  is 
usually  a  coder  layer  nearly  120  m  in  depth  which  is  sta- 
ble, since  it  occurs  within  the  main  halocline.  Below  the 
cooler  layer  there  may  be  an  increase  in  temperature  with 
depth,  an  inversion.  Temperature  inversions  can  be  classi- 
fied into  three  types  according  to  their  method  of  formation 
and  behavior.  The  first  type,  which  is  the  most  common  in 
the  eastern  North  Pacific,  is  a  result  of  the  seasonal  cycle 
of  heating  and  cooling.  The  upper  layers  in  an  extremely 
cold  winter  can  become  colder  than  the  waters  in  the  halo- 
cline below  and  a  positive  temperature  gradient  is  formed. 
With  the  advent  of  the  following  heating  season  the 
temperature  cf  the  surface  waters  increase  causing  a 
(negative)  thermocline  tc  form  above  the  positive  thermal 
gradient.      The    temperature    minimum    thus  generally    is    located 
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Figure   5.       Growth  and   Decay  of   the    Thermccime    at    Ocea 
Station   "P"    (Ficra   Tully   and   Giovando,    1963) 
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just  above  the  halocline  at  120  ra  and  changes  very  little 
daring  the  course  of  the  summer  heating  period.  During  the 
following  winter  season  the  upper  thermocline  decays  and  the 
temperature    inversion  may   be   erased    (Uda,    1963). 

Another  process  which  can  cause  temperature  inversions 
is  the  formation  of  sea-ice  in  the  western  Pacific.  The 
remaining  cold  (-1C) ,  high  salinity  surface  layer  sinks  to  a 
level  of  stability  in  the  existing  halocline  and  is  carried 
eastward  in  the  West  Wind  Drift  where  it  completely  mixes 
with  other  waters  before  reaching  140W.  Depending  en  the 
amount  of  heating  of  the  surface  layers  during  the  transit, 
the  temperature  inversion  may  or  may  not  persist  throughout 
the    year    (Eoden,    1964). 

A  third  process  results  in  a  subsurface  temperature  max- 
imum when  the  warm,  more  saline  Subtropical  Water  intrudes 
under  the  colder  Subarctic  Waters.  The  advecting  water  then 
is  slightly  cccled  and  sinks;  it  remains  warmer  than  the 
surroundings,  however,  and  forms  the  temperature  maximum  at 
depth.  Its  upper  boundary  is  generally  limited  by  the  sta- 
bility criterion   of    the    main   halocline    (Uda,     1963)  . 

As  the  ceding  period  begins  in  late  September,  the  sur- 
face   waters    are      cooled    and   mixed   downward      by   convection    as 
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well  as  by  the  strong  winds  of  the  fall  and  winter  months. 
The  seasonal  thermocline  sinks  and  narrows  until  it  reaches 
the  halocline.  By  late  winder  the  upper  zona  has  reached 
its    minimum   temperature    of  the   year. 

At  depths  below  500  m  temperature  variations  are  small, 
seldom  exceeding  standard  deviations  of  0.04C  (Tabata, 
1960).  Fig.  6  shows  the  typical  temperature  structure  at 
Ocean  Station    "F"  for   both  summer   and    winter    months. 

B.        INTERNAL    WAVE   EFFECTS    CN    THERMAL    STRUCTURE 

In  many  vertical  temperature  structure  analyses,  the 
depths  of  the  isotherms  fluctuate  considerably.  Periodicity 
and  amplitude  cf  the  fluctuations  appear  non-systematic,  but 
result  from  the  superposition  of  a  large  number  of  internal 
wave      trains.  The      limiting   frequencies      range      from      the 

Br unt-Vaisala  tc  inertial;  in  the  Subarctic  Region,  periods 
from  5  minutes  to  at  least  half  a  day  have  been  observed. 
The  vertical  displacement  cf  isotherms  is  of  the  order  cf 
5  m  at  the  top  cf  the  thermocline  and  perhaps  as  large  as 
25  m  at  the  level  of  the  principal  halocline  (Tully  and 
Giovando,     1963) . 
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C.        HORIZONTAL    SURFACE    THEBMAL    STRUCTURE 

Near  the  center  of  the  Gulf  of  Alaska  is  located  an  area 
of  temperature  irinimuni  that  coincides  with  the  dynamic  cen- 
ter cf  the  Alaskan  Gyre.  In  summer  it  roughly  overlies  the 
center  of  the  gyre  as  determined  by  the  salinity  maximum, 
but    in   winter   the  two   are   separated    (Uda,     1963) . 

South  from  the  cold  core  dome  to  approximately  35N  the 
theriral  field  is  oriented  zonally  except  near  the  coastal 
regions.  This  configuration  continues  year  round  with  a 
seascnal  range  cf  7.0C  throughout  the  Transition  Region 
(Uda,  1963).  The  horizontal  thermal  gradient  in  the  north- 
south  direction  is  strongest  in  the  Transition  Zone  where 
gradients    cf    1C/100km  are    common,   as   shewn    in    Fig.    7. 
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Fiaur=  7.   Surface  Mean  Tsmpe  rat.  ure  Structure  For  Seotemb' 
1977  (From  Renner,  1977). 
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IV.       EXPERIMENT 

A.       OBJECTIVE 

For  ocean  monitoring  from  space  to  become  effective,  the 
interior  motions  of  the  oceans  must  be  inferred  from  the 
resulting  signature  of  surface  properties.  The  ability  to 
define  the  subsurface  structure  accurately  by  this  method  is 
becoming    increasingly   utilized  and    understood. 

Legeckis  and  Gordon  (1982)  described  the  ability  of  sat- 
ellite infrared  imagery  tc  locate  sea  surface  temperature 
fronts,  and  the  relation  cf  these  surface  features  to  sub- 
surface temperature  structure.  For  the  3razil  Current  and 
two  warm  cere  eddies,  they  noted  the  mixed  layer  depth  cor- 
related  with    the    surface    temperature   patterns. 

Eernstein,  Ereaker  and  Whrit ner  (1977)  examined  the  eddy 
formations  in  the  California  Current  from  combined  ship,  air 
and  satellite  platforms  and  found  that  the  sea  surface 
temperature  patterns  conformed  to  subsurface  warm  and  cold 
mesoscale      variations.  They   concluded      that      in      unstable 

regions  such  as  the  California  Current,  surface  and  main 
thertcccline  temperature  distributions  have  relationships  to 
each   other.        Thus,       in    such      a    region,      remotely   sensed    sea 
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surface  temperature  patterns  would  be  able  to  interpret  sub- 
surface  circulation. 

From  these  recent  examples  of  satellite-determined  sub- 
surface structure  and  from  other  oceanographic  literature, 
it  appears  that  the  regicns  of  the  worlds  oceans  where 
large  instabilities  occur  (i.e.,  boundary  currents,  coastal 
upwelling  zones,  etc.),  the  subsurface  dynamic  features  ar a 
readily  deduced  from  the  surface  observations.  Areas  of  the 
seas  which  are  quiescent  and  usually  located  in  the  central 
oceans  have  net  been  as  easy  to  describe  accurately.  In 
these  regicns,  the  use  cf  either  dynamical  or  empirical 
methods  could  be  the  only  means  of  forecasting  the  thermal 
structure    cf   the    cceans. 

Numerical  models  are  tased  en  dynamical  principles  and 
will  ultimately  provide  the  basis  for  predicting  oceanic 
processes.  Empirical  methods  have  been  determined  in  the 
past  for  certain  regions  and  for  limited  time  periods,  but 
do    net      appear   to  be      valid  universally.  If   relationships 

could  be  devissd  for  certain  water  masses  or  oceanic 
regimes,  then  the  possibility  exists  of  piecing  together  the 
subsurface  thermal  features  in  the  different  regions  by  uti- 
lizing remote  sensing  to  define  certain  surface 
measurements. 
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Toward  this  goal,  the  first  objective  of  this  thesis  is 
to  determine  what  particular  relationships  exist  between  sea 
surface  temperature  and.  subsurface  thermal  structure  for  a 
relatively  "guiet"  region  in  the  northeast  Pacific  Ocean. 
By  examining  the  climatology,  an  empirical  set  of  equations 
could  then  te  derived  between  sea  surface  temperature,  mixed 
layer  depth,  thermocline  gradient  and  the  depths  of  promi- 
nent   points   within  the   structure. 

The  second  objective  is  to  produce  a  predicted  vertical 
temperature  profile  from  the  best  fit  regression  equations 
for  a  particular  water  region  and  SST  at  a  future  time  and 
different    location.  When  these      profiles    are      compared    to 

actual  temperature-depth  profiles,  accuracy  of  this  statis- 
tical methcd    will  be    determined. 

B.        E&TA    ACQUISITION 

During  the  period  5  to  14  September  1977,  the  U.S.  Naval 
oceancgraphic  survey  vessel  USNS  SILAS  BENT  conducted  ocean- 
ographic  surveys  to  determine  the  thermohaline  and  sound 
velocity  structures  across  the  Califiornia  Current,  the 
Transition  Zcre  and  Subtropical  Water.  For  this  thesis,  the 
only  data  utilized  were  along  a  1155  km  meridional  track 
south  of    Ocean   Weather   Station   "F",    as    displayed   in   Fig.      8. 
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The  data  consisted  of  28  vertical  temperature  profits 
spaced  approximately  42  km  apart;  they  spanned  the  period  9 
to    1  1    September    1977. 

Weather  conditions  during  this  period  were  governed  by  a 
weak  stationary  low  pressure  system  with  mostly  low  cloud 
cover  and  periods  of  light  drizzle  and  showers.  The  wind 
speed  varied  from  5  to  20  knots  from  several  directions 
causing  wave  heights  of  approximately  3  to  6m  along  the 
entire  track    (FNCC,     1977). 

Expendable  bathythermographs  (XET's)  were  used  to  define 
the  thermal  structure.  The  resulting  temperature  traces 
plotted  onboard  ship  were  coded  on  bathythermograph  log 
sheets,  which  were  utilized  as  the  data  points  for  recon- 
struction of  the  thermal  profiles.  Characteristic  of  this 
procedure  is  the  loss  of  the  fine  thermal  structure,  as  only 
prominent  points  en  the  temperature  profile  are  recorded  to 
give   the    general   shape   and   location   of   important    features. 

Another  problem  associated  with  this  method  of  data 
acguisition  is  the  varying  probability  that  the  observer  on 
the  vessel  at  the  time  of  the  X BT  drop  correctly  chose  the 
proper  points  to  reflect  adeguately  the  mixed  layer  depth, 
thermocline   gradient    and    temperature   inversions.         Incorrect 
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Figure   8.      Location    of   XBl   data    Curing   3-11    September    1977 
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recording    cf      the   values    alcng      with   malfunctioning  in      5   to 
10%   cf   the   XETs   can    also    re  a   source   of  potential   error. 

C.       ASSESSMENT    CF   TEMPERATURE    STBUCTURE 

Tc  begin  the  analysis  cf  subsurface  thermal  structure, 
realizing  the  "crude"  nature  of  the  data,  the  temperature- 
depth  points  were  plotted  en  a  computer  graphics  system  and 
the  28  thermal  profiles  were  reconstructed.  From  these  pro- 
files a  vertical  cross- section  of  the  meridional  temperature 
structure  was  drawn,  Fig.  S.  It  was  noted  immediately  that 
numerous  temperature  inversions  existed  along  the  northern 
and  central  reciens  of  the  track,  but  they  ceased  at  approx- 
imately 4  1N.  This  subsurface  border  clearly  defines  the 
Subarctic  Frcnt  which  divides  the  Subarctic  from  Subtrcpic 
Waters;  while  at  the  surface  there  is  no  indication  cf  the 
division.  The  temperature  inversions  varied  in  strength, 
but  generally  shewed  an  increase  in  depth  in  a  southward 
direction. 

The  upper  nixed  layer  displayed  little  variablility  with 
few  transients  found  within  the  structure.  There  were  no 
large  eddies  rct=d  in  the  surface  layer  and  the  mixed  layer 
depth  appeared  relatively  constant  along  the  entire  track. 
Weak      mescscale      features      are   probably      present,         but      the 
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Figure   9.      Meridional  Cross-section   of   Temperature   at    1U5W 
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spatial  resolution  prohibits  their  detection.  The  tharmo- 
cline  was  also  fairly  constant  in  depth  with  the  strongest 
gradients  between  30  to  6C  m .  The  isotherms  in  the  lower 
layers  characteristically  deepened  to  the  south,  with  sev- 
eral areas  of  large  variability  located  near  the  temperature 
inversions.  The  large  warm  intrusion  located  at  U7N  may  be 
caused  by  errors  in  the  data,  as  a  structure  this  size  in 
the  Subarctic  Region  does  not  appear  to  ba  congruent  with 
the  ether  features,  and  is  not  mentioned  in  the  research 
li  teratur e . 

To  assess  further  the  vertical  thermal  structure,  promi- 
nent points,  gradients  and  zones  were  defined  on  a  typical 
temperature  profile  from  the  test  region.  Fig.  10  illus- 
trates the  variables  that  were  used  to  analyze  statistically 
the   thermal    prcfile.      Descriptions   of    these   variables    are: 


•  SST  -  Sea  surface  temperature  was  defined  as  the  sur- 
face temperature  reported  on  the  bathythermograph  logs. 
Since    this    is  the  temperature    recorded   on   the   XBT    trace 


I 

response    time  of  the    XBT   system. 


when    the    XET  Drobe    entered    the   water,      a   slight    dispar- 
ity     f rem    remotely-sensed      SST      may      occur    due      to      the 


•  MLD  -  Mixed  layer  depth  was  visually  determined  as  the 
greatest  depth  where  the  isothermal  or  slightly 
negative  temperature  gradient  in  the  surface  layer 
changed  to  the  stronger  gradient  of  the  underlying 
ther mccline. 

•  BCT  -  Ecttcm  of  seasonal  thermocline  was  visually 
determined  as  the  depth  where  the  seasonal  thermocline 
markedly  changes  into  a  weaker  negative  temperature 
gradient    below. 
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•  TBOT  -  Temperature  at  the  bottom  of  the  seasonal  ther- 
mocline. 

•  DT  -  Difference  in  temperature  between  the  temperature 
at  the  mixed  layer  deptn  and  at  the  bottom  of  seasonal 
thermccline. 

•  DZ  -  Difference  in  depth  between  the  mixed  layer  depth 
and    bcttcm    cf  seasonal  thermocline. 

•  EUZ  -  Depth  of  the  bottom  of  the  main  thermocline  was 
visually  determined  as  the  depth  where  the  thermocline 
changes  to  an  isothermal  or  positive  temperature 
gradient    ^ust  above    the   temperature    inversion   region. 

•  TBUZ  -  Temperature  at  the  bottom  of  the  main  thermo- 
cline . 

•  DTBUZ  -  Difference  in  temperature  between  the 
temperature  at  the  mixed  layer  depth  and  at  the  bottom 
of   main    thermocline. 

•  DZBUZ  -  Difference  in  deDth  between  the  mixed  layer 
depth    and   the  bottom    cf    main   thermocline. 

•  T460  -  Temperature  at  reference  deptn  of  460  m  which 
was  selected  arbitrarily  at  maximum  depth  of  XBT 
reports. 

•  DT5  -  Temperature  change  in  upper  5  m  of  thermocline 
below    the   mixed    layer   depth. 

•  DT15  -  Temperature  channge  in  upper  15  m  of  thermocline 
below   the   mixed    layer   depth. 

•  DT50  -  Temperature  change  in  upper  50  m  of  thermccline 
below   the    mixed    layer  depth. 

•  DT100  -  Temperature  change  in  upper  100  m  cf  thermc- 
cline   kelcw    che    mixed   layer   depth. 

•  DIS    -    Distance    from    northern   most   XBT    drop. 


48 


SST 


DEPTH 


460-- 


DZBUZ 


TEMP 


-  -r-  MLD 


\L 


Figure    10.      Mcdsl  of    Vertical   Temperature    Structure, 
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This  assessment  defines  the  general  structure  primarily 
in    the      upper   layers      above  the      permanent    halocline.  The 

detail  of  the  smaller  structure  such  as  transients  in  the 
mixed  layer  and  multiple  thermoclines  are  not  included 
because  they  are  variable  and  minor  with  respect  to  the 
principal      structure.  The    temperature      inversion      feature 

below  the  main  thermocline  was  not  characterized  as  its 
occurrence  in  the  region  varies  in  location  and  intensity 
during  the  test  period.  In  those  zones  within  the  Subarctic 
Region  where  a  sub-thermccline  temperature  maximum  exists 
permanently  during  the  heating  season,  additional  variables 
should  be    used   tc  define    this   part   of   the   thermal    profile. 

Initially,  the  sea  surface  temperature  gradient  was  used 
as  a  variable  in  the  analysis.  However,  problems  occurred 
in  computing  consistent  gradients  along  the  track  as  the 
calculated  values  were  unrealistic  and  inaccurate  compared 
with  actual  conditions  at  the  time;  this  variable  was 
dropped    frcm   this  study. 

These  varitles  provided  the  statistical  base  from  which 
all  ccrrelaticns  and  regression  equations  were  derived.  The 
values  could  then  distinguish  the  fluctuating  and  transient 
features    from   the  median    thermal      structure.         Tables    I,      II 
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and   III   list      the   values    from   the      analyzed   thermal   profiles 
by   XET  drop   number. 


TABLE    I 
locations    and    Times    of    XBT    Measurements 


XET 

339 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 
351 
352 
353 
354 
355 
356 
357 
358 
359 
360 
361 
362 
363 
364 
365 
366 


CATE/riME(Z) 


090613 

SEP 

77 

47. 

98 

144. 

,  93 

CSC804 

SEP 

77 

47. 

65 

144. 

93 

091003 

SEP 

77 

47. 

20 

144. 

98 

C91200 

SEP 

77 

46. 

82 

144. 

32 

091800 

SEP 

77 

46. 

07 

145. 

05 

092000 

SEP 

77 

45. 

70 

145. 

04 

092205 

SEP 

77 

45. 

33 

145. 

,04 

100000 

SEP 

77 

45. 

02 

145. 

04 

100308 

SEP 

77 

44. 

73 

145. 

,  00 

1005C8 

SEP 

77 

44. 

37 

144. 

98 

100708 

SEP 

77 

44. 

00 

144. 

98 

1CC904 

SEP 

77 

43. 

63 

144. 

98 

101100 

SEP 

77 

43. 

28 

144. 

97 

1G1400 

SEP 

77 

43. 

07 

144. 

95 

101600 

SEP 

77 

42. 

70 

144. 

93 

1C1800 

SEP 

77 

42. 

32 

144. 

95 

102000 

SEP 

77 

41. 

93 

144. 

,93 

11C002 

SEP 

77 

41. 

55 

145. 

00 

110203 

SEP 

77 

41. 

17 

145. 

00 

1104C6 

SEP 

77 

40. 

75 

145. 

00 

110604 

SEP 

77 

40. 

37 

145. 

,  00 

11  1000 

SEP 

77 

39. 

97 

144. 

98 

111200 

SEP 

77 

39. 

60 

145. 

10 

111400 

SEP 

77 

39. 

22 

145. 

25 

111600 

SEP 

77 

38. 

83 

145. 

37 

1  12304 

SEP 

77 

36. 

30 

145. 

47 

120107 

SEP 

77 

37. 

93 

145. 

53 

1203C7 

SEP 

77 

37. 

58 

145. 

67 

LAT  (N)  LONG  (W)  DISTANCE  (km) 


0 

36 

86 

128 

212 

253 

294 

329 

361 

401 

442 

483 

522 

546 

587 

629 

672 

714 

75-7 

803 

846 

890 

931 

973 

1017 

1076 

11  17 

1156 
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TABLE    II 
Values   of   SST,    MLD,    Thermoclir.3   Gradients   and  T460 

XBT      SST  (c)     KID  (l)     DT5(c)    DTJ5(c)     DT50(c)     DT.100  (c)    1460(c) 


339 

13.4 

-33 

1.5 

4.6 

7.1 

6.6 

4.0 

340 

13.6 

-24 

0.3 

1.5 

6.6 

7.4 

4.1 

341 

14.  1 

-25 

2.5 

4.3 

6.0 

6.  2 

5.0 

342 

14.4 

-25 

2.3 

5.6 

7.1 

7.5 

4.2 

343 

15.  1 

-19 

2.4 

4.8 

7.3 

8.  1 

4.6 

344 

15.3 

-32 

0.3 

4.6 

7.7 

7.4 

4.9 

3  45 

15.6 

-29 

3.1 

4.9 

7.3 

7.  9 

6.  2 

3  46 

15.7 

-19 

0.4 

4.7 

6.7 

7.9 

4.7 

347 

16.6 

-27 

0.2 

4.3 

7.9 

9.  1 

4.7 

348 

16.2 

-20 

0.5 

4.7 

7.4 

8.4 

4.8 

349 

17.0 

-21 

0.4 

5.3 

8.0 

8.  5 

4.9 

350 

18.  1 

-20 

0.  2 

6.2 

9.1 

10.3 

4.9 

351 

18.  3 

-30 

3.6 

4.2 

7.8 

10.  2 

5.0 

352 

18.  1 

-20 

0.3 

5.8 

8.4 

8.7 

5.1 

353 

18.4 

-20 

0.7 

4.0 

9.0 

9.9 

5.0 

354 

18.  1 

-20 

0.6 

5.4 

8.6 

9.3 

5.2 

355 

19.3 

-21 

0.6 

5.4 

9.3 

10.  3 

5.  3 

356 

19.7 

-15 

0.2 

6.9 

9.6 

10.8 

5.3 

357 

19.5 

-23 

0.6 

6.4 

9.7 

10.7 

5.3 

358 

20.  1 

-16 

0.3 

2.  1 

9.4 

10.3 

5.6 

359 

20.3 

-22 

0.8 

7.0 

10.0 

10.8 

5.6 

360 

20.  1 

-16 

0.5 

2.6 

8.3 

10.3 

5.9 

361 

20.5 

-30 

2.4 

6.1 

9.3 

10.4 

5.5 

362 

21.5 

-22 

0.3 

5.7 

10.5 

1  1.6 

5.8 

363 

21.4 

-32 

3.8 

6.4 

9.8 

1  1.  4 

5.8 

364 

22.5 

-27 

0.4 

6.2 

9.3 

11.1 

6.5 

365 

22.3 

-22 

0.4 

5.6 

7.6 

9.6 

6.9 

366 

22.7 

-20 

3.3 

5.2 

8.  1 

10.4 

6.8 
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TABLE    III 
Values   of   Thermocline   Tem  pera turss    and   Depths 


III 

EOT (m) 

TECT  (C) 

21  (c) 

11  (1) 

BUZ  (31) 

TBOZ  (c) 

DTBUZ(c) 

DZBUZ  (m) 

339 

-88 

5.4 

7.4 

55 

-95 

5.3 

7.  2 

62 

340 

-30 

6.5 

6.8 

56 

-109 

5.9 

7.4 

85 

34  1 

-70 

7.8 

5.9 

45 

-110 

7.4 

6.  3 

85 

34  2 

-79 

6.6 

7.2 

54 

-123 

6.2 

7.6 

103 

34  3 

-81 

6.6 

8.0 

62 

-138 

6.5 

8.  1 

119 

344 

-81 

7.1 

7.7 

49 

-112 

6.7 

8.  1 

80 

34  5 

-68 

7.8 

7.1 

39 

-112 

7.0 

7.  9 

83 

346 

-57 

8.4 

6.8 

38 

-100 

7.3 

7.9 

81 

347 

-68 

3.6 

7.8 

41 

-114 

7.  3 

9.  1 

87 

348 

-6  5 

8.4 

7.3 

45 

-143 

7.6 

8.  1 

128 

34  9 

-50 

9.1 

7.5 

29 

-109 

3.  1 

8.  5 

88 

350 

-56 

8.8 

9.0 

36 

-123 

7.5 

10.3 

103 

35  1 

-57 

1C.7 

7.0 

27 

-148 

7.  2 

10.  5 

1  18 

35  2 

-73 

9.0 

8.5 

53 

-158 

7.  3 

9.7 

138 

353 

-67 

9.4 

8.8 

47 

-132 

7.7 

10.5 

112 

354 

-63 

9.3 

8.6 

43 

-120 

8.  6 

9.  3 

100 

35  5 

-59 

9.9 

9.1 

38 

-130 

8.  6 

10.  4 

109 

356 

-54 

9.9 

9.5 

39 

-134 

8.5 

10.9 

1  19 

357 

-62 

10.1 

9.1 

39 

-127 

8.4 

10.3 

104 

358 

-50 

11.3 

8.6 

34 

-120 

9.  3 

10.6 

104 

35  9 

-53 

10.8 

9.4 

31 

-95 

9.4 

10.  3 

73 

360 

-78 

10.5 

9.5 

62 

-130 

9.3 

10.7 

164 

36  1 

-99 

9.5 

10.2 

69 

-127 

9.  3 

10.  4 

97 

362 

-68 

11.0 

10.4 

46 

-157 

9.4 

12.0 

135 

363 

-100 

10.2 

1  1  .0 

68 

-155 

9.7 

11.  5 

123 

364 

-91 

12.2 

10.0 

64 

-173 

10.5 

11.7 

144 

365 

-100 

12.5 

9.3 

78 

-196 

10.7 

11.  3 

174 

366 

-100 

12.2 

10.2 

80 

-188 

10.7 

11.7 

168 
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To  examine  the  variations  and  -rends  of  the  individual 
thernal  profile  features  along  the  track,  graphs  of  the 
variables  as  related  to  distance  and  latitude  are  plotted  in 
Appendix  A.  These  diagrams  illustrate  some  of  the  concepts 
that  were  discussed  in  previous  chapters  pertaining  to  the 
variability    in    thermal   structure    within  a    particular   region. 

Ihs  sea  surface  temperature  over  the  entire  track  as 
shewn  in  Fig.  A.1  increases  linearly  to  the  south  where  the 
mean      gradient      was      0.8C/100km.  The     strongest      gradient 

occurred  in  the  vicinity  of  4 4N  and  reached  a  value  of 
2. 3C/100km.  Fig.  11  shews  the  sea  surface  temperature  in 
"-.he  North  Pacific  as  compiled  from  ship  reports  and  clima- 
tology by  Fleet  Numerical  Cceanography  Center  for  10  Septem- 
ber 1977.  In  the  region  of  data  acquisition,  the  meridional 
temperature  gradient  and  surface  isotherm  structure  closely 
compare  tc  the  surface  temperature  versus  distance  plot 
(Fig.      A.1).  Satellite    infrared    imagery    during     this    time 

period  was  obstructed  due  tc  the  thick  cloud  cover  over  the 
area;  therefore,  this  rrethod  of  sea  surface  temperature 
determination  was  not  used  as  was  desired  (Tabata  and 
Kimber,    1579)  . 
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Figure    11.      Pacific    Ocean    Sea   Surface   Temperature   on    10 
September    1977    (Fleet   Numerical  Oceanography 
Center,     1  977)  . 
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Fig-  A.  2  shews  'he  depth  of  the  mixed  layer  along  the 
track.  The  mean  depth  was  -22. 9  ra  with  a  maximum  and  mini- 
mum depth  cf  -32  and  -15  m,  respectively.  The  layer's  gen- 
eral trend  of  a  slight  shallowing  by  about  5  m  in  the 
southward  direction  is  otserved  until  40N  is  reached  where 
the  trend  reverses  and  becomes  deeper.  This  location  of 
change  in  the  mixed  layer  depth  trend  accurately  coincides 
with    the    Subarctic  Front. 

The  fcur  thermocline  gradient  variables  in  Figs.  A. 3 
through  A. 6  displayed  a  measurable  trend  along  the  track. 
The  variability  of  the  temperature  change  in  the  upper  5  m 
of  the  thermccline  (DT5)  was  large  in  the  northern  and 
southern  ends  cf  the  track.  In  the  center,  with  the  excep- 
tion cf  cne  data  point,  a  mean  value  of  0.3C/5m  was 
observed.  The  tcundaries  of  these  three  zones  of  the  BT5 
variable  agreed  with  the  division  of  the  track  into  the 
three  regiens,  i.e.,  Subarctic  Water,  Transitional  Zone  and 
Subtropic  Water.  This  effect  was  noted  in  many  of  the  ther- 
mal structure  variable  plots.  DT15,  DT50  and  DT100  all 
increased  with  decreasing  latitude,  but  the  general  trends 
differed    between     each   variable.  DT15    showed      the   largest 

range     of      variability      with      only        a      slight      increase      in 
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temperature  gradient,  while  DT50  and  DT100  both  increased  at 
a  greater  rat?.  When  plotted  collectively,  as  shown  in  Fig. 
12,      the      relationships    become  apparent.  Generally,       with 

increasing  sea  surface  temperature,  the  negative  temperature 
gradient  of  the  thermccline  increases  slowly  at  a  relatively 
constant   rate. 

Two  ether  depth-dependent  variables  on  the  vertical 
temperature  profile,  besides  the  mixed  layer  depth,  are  the 
bottom  of  the  seasonal  and  main  thermoclines,  BOT  and  EUZ, 
respectively.  Figs.  A. 7  and  A.  1 1  display  these  variables  as 
they  vary  with  latitude.  The  variability  of  the  plots  may 
be  due  in  part  tc  internal  waves  with  typical  amplitudes  of 
15    m    for    the    ECT    and    30    m    for   the    BUZ. 

The  BCT's  general  trend  southward  along  the  track  to 
about  40N  shewed  a  slight  decrease  in  depth  from  -80  to 
-60  m.  Again,  this  location  coincides  with  the  Subarctic 
Front  as  did  the  MLD.  The  BOT  then  increases  sharply  as  the 
trend  reverses   tc  a    depth   of   -90    m. 

For  the  BUZ,  the  trend  is  just  the  opposite  of  the  EOT 
in  the  Subarctic  Water  and  Transitional  Zone.  The  depth  of 
the  bottom  of  the  main  thermocline  increases  slowly  until 
reaching    the    Subarctic  Frcnt,    where    the   maximum    variation    in 
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depth  is  approximately  85  a.  In  Subtropic  Mater,  the  trend 
of  the  BOZ  increases  in  depth  at  a  greater  rate  along  the 
southern    end    cf   the   track. 

When  the  MID,  BOT  and  EOZ  are  plotted  together,  as  shown 
in  Fig.  13,  additional  relationships  appear.  In  the  Subarc- 
tic and  Transition  Regions,  north  of  40N,  the  thickness  cf 
the  upper  layers  above  the  bottom  of  the  seasonal  thermo- 
cline  (BOT)  decreased  in  a  southward  direction.  With  the 
MLD  trend  at  a  relatively  constant  depth,  the  only  thermal 
layer  to  change  in  magnitude  appreciably  is  the  thermocline. 
This  suggests  that  the  seasonal  thermocline,  with  its 
negative  temperature  gradient,  becomes  stronger  as  the  SST 
increases  and  as  the  thermocline  thickness  decreases. 
Indeed,  the  temperature  difference  (DT)  and  thickness  (DZ) 
of  the  seascral  ther mocline  '  displayed  in  Figs.  A. 9  and 
A.  10,   respectively,    agree    with  this    suggestion. 

As  previously  mentioned,  the  dep+h  of  the  bottom  cf  the 
main  thermocline  gradually  deepens  toward  the  south  along 
the  track.  Figs.  A. 13  and  A. 14  show  the  increasing  differ- 
ence in  temperature  and  depth  from  the  mixed  layer  depth, 
DTBUZ  and  DZBOZ,  respectively.  These  tendencies,  in  con- 
junction     with  the      trend      of   the      EOT,         indicate  that      the 
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Figure    13.      Cca-parison   of    MID ,    BCT   and    3UZ    Versus   Latitude, 


60 


negative  temperature  gradient  at  the  lower  part  of  the  main 
thermocline  becomes  weaker  and  thicker  toward  the  south  in 
the    Subarctic   and  Transitional   Regions. 

For  Subtropical  Water,  south  of  40N,  the  BOT  has 
reversed  its  trend  at  the  Subarctic  Front  and  increases  in 
depth.  The  EDZ  continues  to  deepen,  but  at  a  greater  rate 
than  to  the  north.  Temperatures  at  these  two  locations  on 
the  thermal  profile,  TBOT  and  TBUZ,  slowly  increase  linearly 
with  decreasing  latitude  at  nearly  the  same  rate  of 
0.5C/100km.  These        values      suggest      that        the      negative 

temperature  gradient  at  the  lower  part  of  the  main  thermo- 
cline continues  to  weaken  and  thicken,  but  at  a  deeper  depth 
than  in  the  northern  waters.  At  the  upper  layers,  the 
parameters  of  the  seasonal  thermocline,  DT  and  DZ,  increase 
with    increasing    SST. 

D.        STATISTICAL    EROCEDURE 

To  predict  the  subsurface  thermal  structure  from  a  given 
sea  surface  temperature,  a  statistical  approach  was  utilized 
to  determine  the  relationships  between  the  thermal  structure 
variables.  The  data  analysis  was  done  using  the  BIOMED  sta- 
tistical computer  programs  on  the  Naval  Postgraduate 
School's    IEM    3033      computer  system.         Details   of      the    BIOMED 
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computer  programs  used  in  the  analysis  are  provided  in 
Appendix  E.  Easic  statistics  (mean  and  standard  devia- 
tions), scatter  diagrams,  correlation  and  linear  regression 
analysis    were   performed    on  all  variables. 

As  part  cf  the  regression  analysis,  the  correlation 
coefficients  fcr  all  combinations  of  pairs  of  the  sixteen 
variables  were  computed.  These  simple  correlations  are  an 
estimate  cf  the  strength  cf  the  linear  association  between 
two  variables  and  do  not  imply  the  cause  for  the  relation- 
ship  cr    which    variables   affect   the    variation    of   the  other. 

Scatter  diagrams  were  used  in  conjunction  with  the 
regression  analysis  to  assist  in  refining  the  correlations 
between    variables.  By    visually   removing    poor      data    points 

that  detracted  from  trends  within  specific  regions,  the 
analysis  could  be  improved.  After  determining  the  best  fit 
correlations,  least  square  linear  regression  equations  are 
computed  between  variables  and  also  with  sea  surface 
temperature  as  the  independent  variable.  This  enabled  a 
vertical  temperature  profile  to  be  determined  by  calculating 
the  values  fcr  the  variables  with  the  only  input  being  the 
sea    surface   temperature. 
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The  values  cf  the  16  thermal  structure  variables  from 
the  28  temperature  profiles  along  the  track  were  correlated 
for  all  combinations  with  the  results  shown  in  Table  IV. 
Several  strong  relationships  appear  with  the  correlation 
coefficients    cf   a  few   variables    being    greater   than    0.750. 

The  sea  surface  temperature  ha3  the  greatest  number  of 
strong  correlations  which  greatly  enhances  the  idea  of  pre- 
diction cf  the  thermal  prcfile  from  the  3ST.  The  highest 
correlation  coefficient  (0.993)  occurred  when  the  SST  was 
compared      with  the      distance      along      the   track.  This      was 

expected  due  tc  the  relatively  constant  north-south  surface 
temperature  gradient  in  this  region  of  the  Northeast 
Pacific. 

Cther  strong  correlations  with  the  sea  surface 
temperature  were  the  temperatures  at  the  bottom  of  the  sea- 
sonal ana  main  thermocline,  T30T  (0.944)  and  TBUZ  (0.949), 
respectively.  When  plotted  together,  as  shown  in  Fig.  14, 
the  trends  clearly  appear.  The  variables'  similar  relation- 
ships imply  that  the  entire  vertical  temperature  profile 
"leans"  reward  warmer  temperatures  as  the  small  negative 
thermal  gradient  in  the  deeper  layers  increase  with  increas- 
ing   sea    surface    temperature.        The      change   in   temperature   in 
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TABLE    IV 

Correlation  Coefficients    of   Thermal    Variables 


hLD 

SST 

DT5 

DT15 

DT?0 

DT100 

DT 

DZ 

SOT 

TBOT 

3UZ 

TBUZ 

LTBUZ 

DZ5UZ 

DIS 

^60 


MLD 
.221* 
..k8U- 
-.121 
.157 
.235 

.128 
..156 
.■+59 
.296 
..202 
.267 
.217 
•  371* 
.217 
.1  iu 


SST 
.023 
.335 
.770 
.905 
.£35 
.305 
.203 
,91+1+ 
.632 
.9^9 
.960 
.673 
.9^3 
.531 


MEAN        -23.2      18.1 
STD  DEV     5.1        2.b 


DT5 
.080 
.211 
.065 
.026 
.280 
.1+06 
.065 
.12»» 
.035' 
•  068 
.031 
.012 
.'39 

1.2 

1.2 


DT15 
.V35 
•393 
.1+66 
.018 

-.055 
.251 

..098 

•  331* 

.368 

.068 
.365 
.265 

5.0 
1.3 


0T50 

.902  DT100 

.51*3  .867  DT 

-.082  .050  .V25       uZ 

.\2h  .030  -.3^1    -.9^9        BCT 

.632  .322  .688      .150  -.01*1 

-.239  -.u96  -066  -.6*3      .5'1* 

.630  .770  .792      .355  -.2}^ 

.8^2  .961  .391      .20"+  -.11m- 

.2*6  .535  .55!      .578     .WOl 

.738  .373  .ci5     -356  -.25! 

.■+35  ,5f>9  .6"^     .'+05  -.3." 


3.3 
1.2 


9.3        6.5     ^0.3   -72.0 
1.5        '.3     ''*.3     15.' 


TBOT 

•653 

BUZ 

.932 

-.65V 

TBUZ 

.685 

-.635 

.330 

DT3UZ 

.665 

-.9S3 

.661 

.63V 

•  931 

-.633 

•  c63 

•  93^ 

.835 

-.636 

.6  32 

.706 

9.3 

-133.3 

J.1 

9.6 

1.3 

27.  t* 

\.k 

1 .6 

DZBUZ 


.65*+ 

DIS 

.663 

.c^6 

10.2 

530. 9 

31+3.6 

6a 


the  upper   50  and   100  m  cf   the  thermocline   also  displayed 
strong  direct  relax  ions hips  with  the  SST. 

Two  pairs  of  variables  which  exhibited  strong  negative 
correlations,  and  for  similar  reasons,  were  BOT-DZ  (-0.949) 
and  EUZ-DZEUZ  (-0.983).  These  combinations  indicate  that  a 
relationship  exists  between  both  the  SOT  and  BUZ  and  the 
MLD.  This  is  due  to  the  relatively  small  variations  in  the 
mixed  layer  depth  when  compared  to  the  large  fluctuations  of 
the  EOT  and  BCZ.  Therefore,  as  the  bottom  of  the  seasonal 
and  main  t her moc lines  varied,  the  depth  difference  between 
these  variables  and  the  MLD  similarly  changed.  This  could 
provide  an  important  link  in  determining  the  mixed  layer 
depth  if  a  viable  connection  could  be  established  between 
the  EOT  or  BUZ  and  SST.  However,  according  to  the  correla- 
tion coefficients  comparing  the  SST  to  MLD  (0.224),  to  BOT 
(-0.203)  and  to  BUZ  (-0.682),  the  strength  of  these  rela- 
tionships was  weak  as  shewn  in  Figs.  15,  16,  and  17.  Corre- 
lations of  the  ether  variables  with  the  mixed  layer  depth 
revealed  little  as  to  possible  methods  to  determine  this 
important  thermal  structure  variable. 

Due  tc   the  noticeable  differences  in   thermal  structure 
characteristics  among  the  various   regions  encountered  along 
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ths  track,  the  data  ware  next  divided  into  Subarctic,  Tran- 
sitional, Frontal  and  Subtropical  groupings.  These  divi- 
sions were  determined  from  the  previously  discussed  regional 
characteristics  and  thermal  analysis  of  Chapters  II  and  III, 
and  in  retrospect,  they  proved  to  be  an  accurate  reflection 
of    the   boundary    locations   at    the    time    of    data   acquisition. 

The  Subarctic  Water  characteristics  were  -raced  from  the 
northern-most  vertical  temperature  profile  on  the  track  to 
XBT  station  number  345  at  45.33N.  The  Transitional  Zone  was 
observed  from  XET  station  346  to  station  357  at  4  1.17N  fol- 
lowed by  the  Subarctic  Front  between  stations  358  and  359. 
Beginning  at  XBT  station  360  at  (39.97N),  and  continuing 
south  to  the  end  of  the  track,  was  an  ocean  region  contain- 
ing properties  of  the  Subticpic  Hater.  Fig.  18  illustrates 
the  divisions  of  the  track  data  points  into  the  various 
regions. 

To  observe  the  separation  of  the  data  by  groupings  ben- 
tar,  and  to  refine  further  the  correlation  coefficients, 
scatter  diagrams  of  several  of  the  thermal  variables  were 
plotted  using  the  BIOMED  computer  programs.  Since  the  only 
variables  that  correlated  strongly  with  SST  were  other 
temperature-dependent   variables    such   as  TBOT,      DT,      TBUZ    and 
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DTBOZ,  a  closer  inspection  of  the  depth- dependent  variables 
was  desired.  To  ascertain  whether  the  sea  surface 
temperature  can  define  the  depth  of  certain  variables  is 
necessary  in  the  development  of  the  thermal  profile. 
Therefore,  regional  scatxer  diagrams  for  the  three  primary 
depth  dependent  variables  (MLD,  BOT,  BUZ)  and  SST  as  com- 
pared to  the  ether  variables  were  plotted  and  analyzed. 

The  scatxer  diagrams  for  sea  surface  temperature  when 
compared  to  the  ether  variables  are  listed  in  Appendix  C. 
Due  to  the  high  correlation  coefficients  of  these  thermal 
variables  the  trends  of  the  data  points  were  clearly  visi- 
ble. Obvious  linear  relations  occur  between  SST  and  DT50, 
DT100f  DT,  TECT,  TBUZ  and  DTBOZ  through  the  various  regions. 
DT5  and  DT15  also  displayed  an  element  of  linearity  when 
several  of  the  mere  distant  data  points  were  removed.  The 
computer-derived  correlation  coefficients  for  these  two 
variables  were  excessively  low  due  to  -he  remote  points. 
For  the  variable  DT5  shown  in  Fig.  C.2,  the  data  behaved  in 
a  strong  linear  fashion  in  the  Transitional  Region  as  a 
temperature  change  of  0.3C  in  the  upper  5  ra  of  the  thermc- 
cline  was  the  mean.  All  cf  the  thermociine  gradient  vari- 
ables, DT5,  DT15,  DT50  and  DT100,  could  then  be  modeled  by  a 
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linear  equation  related  tc  the  SST,  at  least  in  the  Tran- 
sition  Region. 

Two  scatter  diagrams  which  compare  SST  to  BUZ  and  DZBUZ, 
Figs.  C.  1C  and  C.13,  illustrate  a  perceivable  linear  trend 
between  the  variables  while  the  correlation  coefficients 
reflect  a  moderate  level  of  association,  (-0.637)  and 
(0.678),  respectively.  Again,  with  several  of  the  remote 
data  points  removed,  an  improved  correlation  coefficient 
could  be  computed  which  better  portrays  the  relationship 
with    SST. 

Other  depth-dependent  variables,  for  which  a  visual 
relationship  was  observed  with  SST  in  the  scatter  diagrams, 
were  for  ECT  and  DZ ,  Figs.  C.  6  and  C.9.  Both  exhibited  a 
"V,f  or  wedge  shape  with  data  for  the  Subarctic  and  Subtropi- 
cal Haters  en  either  wing  and  data  for  the  Transitional  Zone 
in   the   center.  If    the    northern   or      southern   regional   data 

were  removed,  then  a  higher  correlation  could  be  obtained 
utilizing  that  from  the  Transitional  Region  alone.  However, 
with  the  visual  trend  in  this  region  closely  following  the 
trend  of  data  from  the  Subarctic  Hater,  data  from  the  north- 
ern end  of  the  track  could  be  used  to  assist  in  the  accurate 
determination      of        relationships      with        the      SST.  The 
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combination  of  data  from  these  two  regions  is  reasonable 
since  the  location  of  the  Subarctic  Front  is  on  the  southern 
boundary  of  the  Transition  Zone,  and  the  water  characteris- 
tics are  siiilar.  Therefore,  the  Subtropic  Water  data 
points  shcuid  fce  removed  to  derive  the  correlations  with 
data   from    only   the  two  northern    regions. 

The  two  data  points  in  the  regional  scatter  diagrams  for 
the  Subarctic  Front  consistently  followed  the  Transitional 
Zone  data  trends.  As  shewn  in  Appendix  C,  for  the  thermal 
variables  EOT,  DZ,  BUZ  and  DZEUZ,  the  Front  provided  the 
sharp  break  between  the  Transitional  and  Subtropical 
Regicns,  tut  in  the  scatter  plots  for  the  remaining  thermal 
variables  the  Subarctic  Front  fell  within  the  trend  of  the 
two   northern   regicns. 

The  scatter  plot  comparing  mixed  layer  depth  to  5ST, 
shewn  in  Fig.  C.1,  initially  appears  not  to  reveal  any  pat- 
tern which  cculd  be  used  to  infer  a  relationship  with  the 
SST.  However,  if  the  Subtropical  Water  data  points  are  not 
included,  a  trend  of  decreasing  depth  is  noticed  and  the 
correlation  coefficient  value  of  0.224  between  these  two 
variables    should   increase    appreciably. 
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The  regional  scatter  diagrams  comparing  mixed  layer 
depth  to  the  ether  variables  did  not  exhibit  any  strong  cor- 
relations and  were  not  included.  However,  several  weak  pat- 
terns were  observed  which  indicated  a  trend  of  the  MLD 
slowly  decreasing  in  depth  in  a  southward  direction  for  the 
northern  twe  regions.  All  scatter  plots  displayed  large 
variability  ameng  the  data  points  which  agrees  with  the  lew 
correlation   coefficients    for   the    variable. 

Appendices  D  and  E  contain  the  regional  scatter  diagrams 
for  ECT  and  EUZ,  respectively.  These  plots  show  a  greater 
linear  dependence  among  the  variables,  and  therefore  higher 
correlation  coefficients,  than  was  shown  by  the  mixed  layer 
depth.  In  most  cases  when  the  Subtropical  Water  data  are 
removed,  leaving  only  data  from  the  Subarctic,  Transition 
and   Frontal   regions,    a   stronger    relationship   was   observed. 

Two  noteworthy  scatter  diagrams  which  had  strong 
negative  correlations  were  BOT-DZ  and  BUZ-DZBUZ  as  shewn  in 
Figs.  19  and  20.  As  previously  noted,  these  two  variables 
with  high  correlation  coefficients  could  provide  the  neces- 
sary link  to  determine  the  mixed  layer  depth.  Since  DZ  and 
DZBUZ  are  the  differences  in  depth  from  the  MLD  to  the  EOT 
and  EOZ,  respectively,  a  pair  of  simple  equations  define  the 
relationship. 
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MLD    =    EOT    ♦     DZ 
MLD    =    EDZ    ♦     CZBUZ 
note:    MLE,    EOT    and    BUZ  are    negative   values    (depth) 

The  only  weak  "link"  in  the  two  equations  is  the  deter- 
mination of  the  EOT  and  BCZ  directly  from  SST.  To  improve 
the  correlation  coefficients  for  these  two  variables  as  sug- 
gested by  the  regional  scatter  diagrams,  data  from  the  Sub- 
tropical Water  and  several  distant  data  points  which 
detracted  from  the  established  trends  were  removed.  In  sev- 
eral cf  the  clots,  primarily  for  the  temperature-dependent 
variables,  where  the  trend  was  actually  stronger  and  more 
accurate  with  data  from  the  Subtropical  Region  included,  the 
correlation  coefficients  were  improved.  The  correlation 
coefficients  for  the  best-fit  trends  are  listed  in  Table  V. 
Indications  are  made  to  shew  which  coefficients  increased  or 
decreased    frcu   the   values    presented   in   Table    IV. 

An  example  cf  a  correlation  coefficient  that  decreased 
due  to  the  removal  of  data  points  and  yet  provided  a 
stronger  treed  was  the  KID  versus  DT5  combination  which 
changed    from    -0.484    to   -0.197.       In    this   case,       the      majority 
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Figure    19.      Regional    Scatter   Diagram   -   DZ    Versus   30T 
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Figure  20.   Regional  Scatter  Diagram  -  DZBUZ  Versus  BUZ 
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TABLE    V 
B€St-Fit   Correlation   Coefficients   of  Thermal    Variables 


KLD 

MLD 

SST 

.535' 

SST 

DT5 

0 
-.197 

.271' 

DT5 

DT15 

.197 

.677 

0 
-.063 

d:i5 

DT50 

M2 

.930' 

-M5 

0 
.+60 

DT50 

DT100 

A95 

.905 

-.239 

.335° 

.926 

DT100 

DT 

M9 

-.^97 

.51/ 

.952' 

.8W° 

DT 

DZ 

-.1  56 

-.6V5 

0 
.030 

-.223' 

-.M*1 

-.559  ■ 

-.200 

DZ 

30? 

1 

.557 

.776 

9 
-.255 

.256 

.59* 

.725' 

.4-10 

-.9»*9 

isoT 

• 

T30T 

.Wo' 

.9M+ 

-.199' 

.225 

.77/ 

.839 

0 
.555 

_    779 

1 

.373 

"   ■    1  "• 

3UZ 

-.307' 

-.790* 

.212 

.223' 

-.239 

-.W96  ■ 

-.  5o6 

-.6*3 

.5i»* 

-.653 

BUZ 

T3UZ 

1 

.622 

.9*+9 

-.379' 

.3* 

.71  f 

9 
.fc73 

-.62:+ 

.331' 

0 

.578 

-  651* 

TBUZ 

DT3UZ 

1 

.V36 

.960 

-.30"+' 

•  357 

1 

.935 

1 
.977 

0 

.359 

-.503 

.650 

.C-+1 

-.63c 

.693° 

DTSUZ 

DZ3DZ 

.531 

/ 

.766 

.021 

/ 

.120 

.221 

.293° 

.292 

.237 

0 
.171 

.132" 

-.933 

.071° 

0 
oW 

DIS 

/ 

.55* 

.993 

,hQ? 

.6W 

.917 

/ 

.922 

.355 

1 

-.633 

.775' 

.931 

-.665 

.963 

.93'+ 

.cc* 

DIS 

TWO 

.2^0 

.331 

.057 

.276 

O-l 

.599 

,6Mt 

-.  5S3 

.598 

.535 

-.686 

.882 

.7C6 

•  6  3 

.8W 

not€ :    (M     indicates    increased  coefficients 
(°)     indicates   decreased  coefficients 
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of  data  pciEts  were  grouped  about  the  DT5  value  of  9.4C  per 
5  m  and  a  MLC  value  of  -23  m.  The  attempt  to  capture  the 
trend  of  the  Transition  Fegion  was  made  by  excluding  the 
distant  scattered  points;  however,  this  detracted  from  the 
overall   trer.d    and  reduced    the   correlation    coefficient. 

A  large  number  of  correlation  coefficients  improved  sig- 
nificantly by  this  method  cf  data  exclusion.  Most  notice- 
able is  the  variable  SST,  which,  (except  for  the  combination 
with  DT5)  increased  all  ether  coefficients  to  greater  than 
0.645  with  six  ircreasing  above  0.900.  The  important  corre- 
lations of  SST  with  BOT  and  BUZ  improved  from  -0.203  and 
-0.682  to  0.776  and  -0.790,  respectively.  The  correlation 
coefficient    cf    SST   versus    EOT    even    reversed   sign. 

This  increased  strength  in  the  various  correlations 
between  variables  provided  the  confidence  to  infer  that  usa- 
ble relationships  exist  in  the  vertical  thermal  profile. 
Since  the  variable  SST  had  the  greatest  number  of  strong 
correlations  with  the  other  variables,  the  statistical 
method  of  deriving  linear  regression  equations  using  only 
the    input    cf    SST    appeared   practical. 
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E.        REGRESSION    EQUATION    ANALYSIS 

From  the  rest  fit  correlation  coefficients  of  the  vari- 
able SST  in  combination  with  the  other  variables,  the  fol- 
lowing  linear    regression    equations    were  derived: 

HLD    =     (1.2C01)     SSI    -    43.300 

ET5    =     (0.03629)     SST    -    0.21573 


DT15    =     (0-16058)     SST    ♦    2.2271 


DT50    =     (0.48462)     SST    -    0.14105 


DT100    =     (0.48506)      SST    +    0.52999 


BCT    =     (3.9611)     SST    -    133.0 


TBOT     =     (0.61762)      SST    -     1.9317 


DT    =     (0.41729)     SST    +    0.94677 


DZ    =     (-2.7611)     SST    +     89.782 


EUZ    =     (-7.5247)     SST    -    0.00793 


TEUZ    =     (0.48510)      SST    -     0.66001 


DTEUZ    =     (0.55345)     SST    -    0.42130 


DZEUZ    =    (7.7697)      SST    -27.820 
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To  construct  a  complete  vertical  temperature  profile, 
since  the  variables  previously  mentioned  only  defined  the 
upper  domain  cf  the  therial  structure,  a  new  variable  was 
needed  from  lower  depths.  The  temperature  at  460  in  was  cho- 
sen because  it  was  the  lowest  datum  point  available  from  the 
XBT  data  provided  along  the  track.  This  type  of  variable 
was  ncc  particularly  desired  as  it  relates  a  temperature  tc 
a  specific  depth  rather  than  to  a  feature  en  the  thermal 
profile.  A  better  variable  would  be  at  a  reference  decih 
where  the  t eipeiat ure  gradient  becomes  isc-hermal  in  the 
deeper  regions  of  *he  ocean.  However,  this  was  not  availa- 
ble from  the  United  depth  of  the  data  and  the  T460  variable 
was  established.  The  same  methods  of  analysis  were  employed 
on  this  variable  as  were  conducted  on  the  or  hers.  When  ccm- 
oared  to  the  sea  surface  temperature,  a  correlation  coeffi- 
cient  of  0.321  was  computed  along  with  a  strong  linear 
trend  in  the  regional  scatter  diagram.  Therefore,  a  linear 
regression  equation  could  be  used  to  determine  the  bottom 
point  of  the  predicted  temperature  profiles  at  460  meters 
from  the  variable  SST.  The  linear  regression  equation  for 
T460    was    computed  as: 

T460    =     (0.21834)     SST    +     1.3109 
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Two  ether  important  relationships  which  were  derived  by 
regression  analysis,  but  they  do  net  use  SST  as  the  indepen- 
dent   variable: 

DZ  =     (-0.8539U)     BCT    -    12.693 

DZBUZ    =    (-1.0359)     3UZ    -    27.392 

These  equations  use  the  computed  values  of  BOT  and  3U2 
as  determined  from  the  SST  to  calculate  the  thickness  of  the 
seasonal  and  main  thermcclines .  The  magnitude  of  these 
thicknesses  (IZ ,  DZEUZ)  when  added  to  the  BOT  or  BUZ, 
respectively,    can  then   define   the   mixed  layer    depth. 

Due  to  the  relatively  low  correlation  coefficient 
between  SST  and  MLD  for  the  SST-derived  MLD  equation 
(0.535),  possibly  a  combination  of  equations  which  define 
3LD  with  higher  correlations  could  improve  the  accuracy  in 
determining  the  mixed  layer  depth.  The  regression  equation 
for  the  SST-derived  BUZ,  with  a  correlation  coefficient  of 
-0.790  between  the  two  variables,  along  with  the  3UZ-derived 
DZBUZ  mentioned  above  with  a  correlation  coefficient  of 
-0.983,  seemingly  provide  the  strongest  combination. 
Therefore,  the  group  of  equations  used  for  the  mixed  layer 
depth   determination    were: 


BDZ   =     (-7.5247)     SST    -    0.30793 
DZEDZ    =     (-1.0359)     BUZ    -    27.392 
MLD   ■    BUZ    +    DZBUZ 

Three  ether  combinations  of  regression  equations  could 
have  been  used  to  determine  the  mixed  layer  depth,  but  were 
not  chosen  because  of  their  lower  correlation  coefficients 
between    variables.  The    equations   and      correlation    coeffi- 

cients   (r)    are    listed   belc*: 

BOT    =     (3.S611)     SST    -     133.08  r    =    0.776 

DZ    =     (-2.7611)     SST    ♦    89.732  r    =    -0.645 

HLD    =    BCT    +    EZ 


BOT    =    (3.S611)     SST    -     133.08  r    =    0.776 

DZ    =     (-0.65394)     EOT    -    12.693  r    =    -0.949 

MLD    =    BCT    ♦    DZ 
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BUZ    =     (-7.5247)     SST    -   0.00793  r    =    -0.790 

DZBUZ    =     (7.7697)     SST    -    27.820  r    =    0.766 

MLD    =    BUZ    ♦    DZBUZ 

With  the  MID  determined,  the  tempera ture  at  this  depth 
was  needed  tc  provide  the  fundamental  temperature  from  which 
the  tharmccline  could  be  computad.  A  variable  TMLD  was 
added  to  represent  the  temperature  at  the  mixed  layer  depth 
and  was  defined  as  the  sum  of  TEUZ  and  DTBUZ.  Since  both 
temperature-dependent  variables  had  very  high  correlations 
with  the  SST,  0.949  and  0.960  respectively,  strong  confi- 
dence  could    be   placed   in    the    resulting   variable,    TMLD. 

An  additional  feature  in  calculating  TMLD  by  this  method 
was  to  allow  for  a  slight  negative  temperature  gradient  tc 
occur  in  the  ttixsd  layer.  Many  ocean  thermal  structure  mod- 
els assume  that  this  upper  layer  will  be  isothermal  when 
actually  a  realistic  profile  often  includes  a  temperature 
gradient. 

The  four  thermal  gradient  variables  DT5,  DT15,  DT50  and 
DT100,  as  determined  from  the  SST,  were  subtracted  from  the 
temperature  and  depth  of  the  mixed  layer  to  form  the  thermo- 
cline.        This    structure,      along    with   the    SST    and   mixed    layer 
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depth  values,  provided  the  upper  portion  of  the  vertical 
temperature  profile. 

In  the  region  below  the  seasonal  thermcclir.e  the 
temperature-dependent  variables,  TEOT  and  IBUZ,  and  depth- 
dependent  variables,  BOT  and  BOZ,  were  computed  directly 
from  the  input  cf  SST  and  applied  to  the  temperature  pro- 
file. Several  times  during  the  plotting  of  these  points  in 
the  testirg  and  ccmpariscr  phase  of  the  thesis,  the  location 
of  the  DT100  variable  would  occur  below  both  the  BOT  and  BUZ 
or  between  them.  In  these  cases  it  was  simple  to  adjust  the 
sequerce  cf  points  to  allow  the  proper  construction  of  the 
temperature  profile. 

Finally,  in  the  generation  cf  the  temperature  profiles, 
the  variable  TU6C  was  used  to  provide  the  lower  base  of  the 
structure.  Since  TU60  was  an  artificially-generated  vari- 
able, defined  at  a  constant  depth,  the  characteristic  small 
temperature  variations  at  that  depth  should  still  be  accu- 
rately reflected. 

A  total  cf  nine  points  were  used  to  define  the  entire 
vertical  temperature  profile.  These  were  plotted  on  a 
TSKTHCNIX  618  graphic  display  system  and  connected  with  a 
rational  spline   method  that  allows  some   curvature  between 
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points.         This   was   done   to     portray   the  temperature   profiles 
of   the  regicn   mere   realistically. 

F.        THERMAL    FFCFILE    TEST    AND    COMPARISON   ANALYSIS 
1.       TJST    J 

The  first  type  of  comparison  test  was  conducted  with 
the  original  -rack  of  data  used  in  the  generation  of  the 
predicted  thermal  profiles.  The  XET  temperature-depth  pro- 
files of  the  Subarctic  Water  and  Transition  Zone  were  com- 
pared with  the  computed  temperature  profiles  which  utilized 
the  sea  surface  temperature  of  the  actual  XET  trace  as  the 
only  input.  Fig.  21  shows  the  comparison  between  the 
observed  and  predicted  thermal  structures  for  the  northern- 
most XBT  station  on  the  track.  Appendix  F  lists  the  remain- 
ing twenty  profile  comparisons  for  the  northern  regions 
excluding    the    Subtropic    Water. 

The  comparisons  generally  illustrated  a  strong  simi- 
larity in  the  upper  layers.  The  mixed  layer  depth  along 
with  the  slight  negative  temperature  gradient  in  the  layer 
were  realistically  predicted  for  a  majority  of  cases;  inter- 
nal wave  motions  are  assumed  responsible  for  the  largest 
variation  from  the  predicted  profiles.  The  upper  thermo- 
cline     region        closely      approximated      the        strong     thermal 
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gradient    measured  by      the    XBTs.         In   the      lower   thermccline, 

greater  variability  occurred,  assumed  due  to  internal  waves 
and  therial  advection,  resulting  in  poor  comparisons  as 
shown  in  Figs.  F.2,  F. 1 1  and  F.  12  for  XBT  stations  341,  350 
and    351. 

The  predicted  lower  regions  were  not  as  precise  when 
compared  to  the  observed  profiles  usually  due  to  the  exis- 
tence of  a  significant  temperature  inversion.  Lack  of  a 
model  for  this  feature,  which  was  not  attempted,  clearly 
affected  the  profiles  at  depths  greater  than  130  m.  Several 
predicted  profiles,  however,  did  attempt  to  reproduce  this 
structure  with  the  placement  of  the  DT100  variable  below  the 
BUZ.  This  was  shown  in  XET  sta  tions  343-346  and  348  (Figs. 
F.  4    -   F.7    and    F.9)  . 

Between  the  bottom  of  the  main  thermocline  and  460 
m,  several  predicted  profiles  completely  failed  to  approxi- 
mate the  observed  structure.  This  was  caused  by  the  lack  of 
variables  in  the  lower  layers  as  well  as  the  computer-gener- 
ated curvature  between  data  points.  When  the  variables 
DT100  and  BUZ  were  very  close  together,  as  shown  for  XET 
stations  348  and  349  (Figs.  F.9  and  F.10),  the  profile  wculd 
bulge   with   a    positive   or    negative   gradient. 
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Utilizing  only  nine  points  to  define  the  predicted 
thermal  profile  obviously  smooths  the  fine  structure  within 
the  thermocline  and  mixed  layer.  However,  due  to  the  small 
size  and  transient  nature  cf  these  features,  their  removal 
would  not  severely  limit  the  usefulness  of  the  predicted 
structures . 

In  this  initial  test,  since  the  predicted 
temperature  structures  were  generated  and  compared  with  the 
same  XBT  data  from  the  meridional  track,  the  comparisons 
should  be  relatively  close.  The  predicted  profiles  would 
actually  represent  the  average  of  all  the  observed  thermal 
structures  as  the  variables  reflect  the  trends  along  the 
track. 

2.   TEST  2 

A  more  stringent  and  realistic  test  was  attempted 
that  defines  the  spatial  limitations  of  this  thermal  pre- 
diction method.  Test  2  consists  of  comparing  observed  and 
predicted  temperature  profiles  which  occurred  during  the 
same  time  period  as  the  data  generation,  but  spatially 
removed.  A  possible  indication  as  to  a  limit  in  the  effec- 
tive distance  would  be  the  observation  that  the  predicted 
profile  begins  to  diverge  greatly  from  the  actual 
temperature    profile. 
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The  test  XBT  data  were  obtained  from  Fleet  Numerical 
Oceanography  Canter  by  a  computer  search  for  all  available 
XBT  drops  which  cccurred  in  September  from  1965  to  1980  in  a 
ten  degree  wide  strip  of  longitude  which  surrounds  the  orig- 
inal track.  Inexplicably,  the  data  received  were  limited  to 
the  area  south  of  44N,  excluding  the  Subarctic  Water.  Lack 
of  tine  prevented  an  attempt  tc  acquire  these  data  whicn 
would  be  useful,  but  not  entirely  necessary,  for  analysis. 
The  temperature-depth  data  were  then  selected  by  noting  the 
positions    relative      tc   the     original   data.  To   provide      an 

effective  test,  a  variety  of  locations  in  the  Transitional 
Zone  were  chosen.  However,  these  positions  were  constrained 
by  the  availability  of  data  in  the  region  and  appeared  as  a 
zonal  belt  perpendicular  tc  the  original  track.  Fig.  22 
shows   the    location  of   original  and    test   XBT   drops. 

The  temperature-depth  pcir.ts  proviied  by  FNOC  were 
then  plotted  and  connected  on  the  computer  graphical  display 
with  the  same  amount  of  curvature  applied  tc  the  profile  as 
before.  Ecth  observed  and  predicted  profiles  were  plotted 
together    again    for  analysis  and    comparison    purposes. 

Five  vertical  temperature  profiles  from  the  10-11 
September    period     were   compared    with   the      predicted   profiles 
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as  generated  from  the  sea  surface  temperature  of  the 
observed    profiles.  The    test    produced   varied     results    with 

several  features  of  the  computed  temperature  profiles  com- 
paring favorably  with  the  actual  thermal  structure.  The 
largest,  differences  were  noted  for  the  two  profiles  located 
the  greatest  distance  frcm  the  meridional  track:  at  145*. 
Figs.  23  and  2U  illustrate  the  comparisons  at  these  loca- 
tions, approximately  400  km  east  and  west  from  the  original 
track  of  data.  The  mixed  layer  depth  and  upper  thermocline 
regions  cf  the  predicted  profile  closely  resembled  the 
observed  profile,  but  both  structures  displayed  differences 
in   the   lower    regions. 

The  eastern  most  observed  profile  had  a  thin,  strong 
thermcclire  which  gradually  weakened  with  depth.  This  con- 
trasted markedly  with  the  strong  and  sharply-defined  thermo- 
cline of  the  predicted  profile.  The  regions  below  150  m 
displayed    similar  temperature    gradients   for    both   profiles. 

Fcr  the  western  meet  test  profile,  the  only  region 
to  differ  significantly  frcm  the  predicted  structure  was  the 
lower  domain  below  130  m.  The  location  of  the  two  close 
variables,  BDZ  and  DT100,  along  with  the  computer -gen erated 
curvature  forced  the  predicted  profile  to  diverge  from  the 
weak    temperature    inversions. 
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Figure   23. 


Observed    and    Predicted    Thermal   Profiles   a.- 
41. 5N,     140. 3W    For    101200Z    SS?    77. 
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Figure   24.      Observed    and    Predicted   Thermal   Profiles  a' 
42. 3N,    150.0M    For    110059Z    SSF    777 
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The  three  other  test  cases,  shown  in  Figs.  25,  26 
and  27,  which  were  located  closer  to  the  original  track,  all 
displayed  strong  similarities  between  the  observed  and 
predicted  profiles.  Except  for  the  fine  structure  and  pos- 
sible effects  of  internal  waves  en  the  depr.h  of  the  thermo- 
cline,        the   computed      profiles      were      very    accurate.  The 

unusual  shape  of  the  observed  tempera-cure  profile  at  the 
bottom  of  the  seasonal  theimocline,  in  Fig.  27,  and  close 
approximation  by  the  predicted  profile,  illustrates  the 
strong  ability  of  this  statistical  method  to  produce  thermal 
features    unique    to   a    certain   region. 

To  test  the  spatial  limitations  of  this  empirical 
method  further,  a  vertical  temperature  profile  taken  at 
Ocean  Weather  Station  "P"  during  the  data-generation  time 
pericd  was  examined.  Since  the  location  of  the  northern 
most  X3T  drop  en  the  track  was  two  degrees  of  latitude  south 
of  OWS-P,  the  results  should  indicate  the  reliability  of  the 
method  outside  the  test  region.  Fig.  28  displays  the  com- 
parison between  the  observed  and  predicted  profiles  for 
OWS-E. 

Frcm  the  surface  tc  the  middle  of  the  thermocline, 
the    two      profiles  were   very     similar,      with   the      mixed    layer 
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Figure    25.      Observed    and    Fredicted    Thermal    Profiles   a* 
U1.9N,    146.7W    for    101200Z    SEP    77. 
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Figure    26. 


Observed    and    Eredicted   Thermal    Profiles   a- 
42. 3N,    144. OW    For     102115Z    S2F    77. 
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Figure    27. 


Otserved    and    Predicted   Thermal    Profiles   aJ 
U1.1N,     143. 2W    For    110000Z    SEP    77. 
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Figure   28.      Cts«=rved    and    Eredict€d   Thermal    Profiles   at    CWS-? 
fcr     1103182    SEE    77. 
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depth  accurately  depicted.  However,  the  lower  segment  of 
the  two  profiles  differed  by  2.0C.  With  the  actual  vertical 
temperature  profile  being  considerably  cooler  than 
predicted,  a  different  regime  must  have  been  encountered. 
Tabata  (1961)  suggests  ficm  his  observations  in  the  region 
that  horizontal  transport  cf  colder  water  from  the  vicinity 
of  the  Alaskan  Gyre,  lying  northwest  of  the  station,  may 
cause  a  decrease  in  temperature  during  the  summer  months. 
3.       TEST    3 

The  third  test  included  bo-h  spatial  and  time  dif- 
ferences in  an  attempt  tc  forecast  the  vertical  thermal 
structure  up  to  three  weeks  in  the  future.  The  only  availa- 
ble XET  test  da*a  in  the  Transitional  Region  were  located 
within  three  degrees  longitude  to  the  west  cf  the  original 
track  as  secwb  ir.  Fig.  29.  Vertical  temperature  profiles 
were  chosen  at  weekly  intervals  and  compared  with  the 
predicted    structures. 

At  the  er.d  cf  the  first  week,  the  observed  thermal 
profile,  as  shewn  in  Fig.  30,  displays  a  deeper  mixed  layer 
to  a  depth  of  4C  m  cverlying  an  extremely  strong  negative 
temperature  gradient.  This  feature  was  net  represented  well 
by    the    predicted    profile    as  the    computed    MLD    and   thermocline 
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Figure  29. 


Lo  c  a t  io  n 

Data. 


of  Original  (•)  and  Test  3  (X)  XBT 
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variables      were      shallower     in   depth.  The      shape      of      the 

observed  mixed  layer  infers  strong  wind  and  convective  mix- 
ing, so  surface  weather  charts  for  the  period  of  11-18  Sep- 
tember 1977  in  the  Northeast  Pacific  were  checked  fcr  ths 
occurrence  of  a  significant  atmospheric  event.  On  17  Sep- 
tember 1977,  a  cold  front  associated  with  a  low  pressure 
system  in  the  Gulf  of  Alaska  had  quickly  passed  through  the 
test  region  with  rain  and  35  knot  winds.  The  weather  during 
the  period  cf  data  acquisition,  8-11  September  1977,  as 
described  befcre,  was  mcstly  cloudy  with  periods  of  light 
winds  and  scattered  showers  (FNOC,  1977).  These  contrasting 
atmospheric  situations  suggest  a  major  limitation  in  the 
statistical  model,  as  weather  effects  on  the  upper  vertical 
temperature  structure  must  be  prcperly  accounted  for  in  the 
predicted    profile. 

The  lower  half  of  the  thermccline  on  the  predicted 
temperature  profile  closely  followed  the  observed  struc^urs 
down  to  170  m.  Except  for  the  weak  temperature  inversion  at 
200  i,  which  was  not  included  in  the  empirical  model,  the 
lower    domain    tehaved    well. 

The  next  test  profile  was  taker  on  24  September 
1977,        thirteen    days      after      the      data   generation      for      the 
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Figure   30. 


Observed    and    Predicted    Thermal    Profiles   fo' 
182030Z    SEP    77   at    43 . 7N,  146. 1 W . 
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predicted  structure.  Location  of  this  test  XBT  was  on  the 
southern  boundary  of  the  Transitional  Zone  approximately  200 
km  west  of  the  original  track.  The  predicted  temperature 
profile  in  Fig.  31  clearly  resembled  the  observed  vertical 
theriral  structure,  but  with  slightly  cooler  temperatures. 
Due  to  the  strong  winds  from  the  storm  a  week  earlier,  the 
observed  MID  remained  deep.  However,  a  transient  at  30  m 
depth  had  formed  with  a  0.5C  increase  in  sea  surface 
temperature  as  the  weather  in  the  region  during  the  second 
week  was  dominated  by  a  strong  high  pressure  system.  The 
observed  profile's  thermocline- a  Isc  deepen  2d  by  about  10  m, 
which  may  be  due  to  internal  waves,  or  it  may  have  been 
caused  by  wan  advection  from  the  south.  Immediately  below 
the  thermocline,  the  two  profiles  differed  slightly  because 
of      the      fine      structure      in      the      region.  The      predicted 

temperature  at  460  m  depth  remained  accurate  when  compared 
to  the  observed  temperature  throughout  all  the  temporal 
tests. 

The  final  test  forecast  was  conducted  for  an  X3T 
location  209  kir  to  the  west  of  the  north-south  track  on  30 
September  1977,  nearly  two  weeks  after  the  basic  data  were 
acquired.      Another  cold    front   associated   with   a    low    pressure 
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Figure    31.      Otserved    and    Predicted   Thermal    Profiles 
2<4C0CGZ    SEP    77    at    41. 2N,    147. 5tf. 
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system  moved  through  the  area  en  25-26  September  1977  and 
was    immediately      followed    by     a    high      pressure    ridge.  The 

storm,  with  iraximum  winds  of  40  knots,  maintained  the  40  m 
thick  isothercral  mixed  layer  as  shewn  in  Fig.  32-  The  wind 
mixing  had  increased  and  deepened  the  negative  temperature 
gradient  cf  the  thermocline  to  a  level  approximately  10  id 
below  the  predicted  profile.  In  the  lower  domain  of  the 
thermal  structure  a  large  difference  was  noted  due  to  the 
slight  temperature  inversion  of  the  observed  profile  and  to 
the  close  location  of  the  two  variables,  3UZ  and  DT100,  at 
130  n:  on  the  predicted  profile.  The  profiles  coincided  at 
the    460    m    depth   as  determined   by    the   variaole   T460. 

The  accuracy  of  the  predicted  profile  after  three 
weeks  from  the  time  of  data  acquisition  for  the  statistical 
analysis  shows  an  overall  degradation.  However,  considering 
that  two  ston  systems  passed  through  the  region  and  the 
summer  heating  season  was  approaching  its  end,  the  fore- 
casted  profile    was  not   excessively    poor. 

G.        REVIEW    CF    EREDICTED    TEERMAL    PROFILE    RESULTS 

The  underlying  requirement  of  the  statistical  approach 
in  predicting  the  vertical  temperature  profile  is  to  acquire 
accurate    data    that    reflects  the      mean    thermal    conditions    for 
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a  particular  season  and  region.  Besides  the  acquisition 
techniques,  the  importance  in  precisely  defining  the 
regicnai      boundaries    is      essential.  Both  temperature      and 

salinity  shculd  be  used  tc  establish  the  limits  of  usable 
data  when  applying  correlation  and  regression  analyses. 
Results  of  the  predicted  thermal  profiles  located  outside 
the  data  generation  area,  as  defined  by  a  certain  region, 
displayed  a  marked  degradation  in  the  lower  thermociine  and 
lower  domain  regions.  The  predicted  profiles  located  near 
the  area  generally  portrayed  the  observed  temperature  pro- 
files realistically  within  the  spatial  limits  of  455  km 
about   the    criqinal  data. 

The  results  of  the  temporal  accuracy  of  the  predicted 
profiles  indicate  a  favorable  ability  to  forecast  the  ther- 
mal structure  with  certain  limitations.  The  main  limiting 
factor  was  the  weather  which  must  be  included  in  the  analy- 
sis. Clearly,  when  high  winds  pass  through  a  region  after 
the  vertical  temperature  data  is  acquired,  the  upper  thermal 
structure  will  be  altered.  Two  regions  of  the  thermal  pro- 
file that  will  be  affected  are  the  depth  of  the  mixed  layer 
and  strength  of  the  thermociine.  For  short  periods  of  about 
one      week,        this  upper      domain      in      the    water      column      will 
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Figure    32. 
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primarily  reflect  changes  due  to  "he  wind  and  free  correc- 
tive mixing.  Fcr  time  scales  cf  greater  than  a  week,  an 
occasional  thermal  intrusion  from  the  Subarctic  or  Subtrcpic 
Region  may  te  an  additional  factor  than  could  alter  the 
therial  structure.  Such  advective  processes  would  be  diffi- 
cult to  forecast  by  the  statistical  method  in  regions  where 
large  meanders  and  eddies  occur,  but  they  should  net  greatly 
affect  a  relatively  quiescent  area  such  as  in  the  Tran- 
sitional   Zone. 

An  exauple  cf  another  problem  encountered  oy  this  method 
of  thermal  structure  prediction  from  the  sea  surface 
temperature  is  the  process  of  heating  at  the  ocean's  sur- 
face, called  the  "afternoon  effect."  Fig.  33  illustrates 
the  difference  in  the  predicted  profile's  structure  with  the 
observed  "afternoon  effect."  Since  the  sea  surface 
temperature  determines  the  thermal  profile  in  this  statisti- 
cal model,  an  error  is  imposed  in  the  subsurface  structure, 
throughout      the    entire      depth.  With      the    surface      heating 

removed,  the  predicted  profile  wculd  closely  resemble  the 
observed  profile.  The  "afternoon  effect"  would  greatly 
hinder  the  ability  of  remcte  sensing  to  determine  the  sub- 
surface  thermal    structure    fcy    this    method. 
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Ihe  method  of  determining  mixed  layer  depth  was  reviewed 
to  ensure  that  the  best  possible  value  was  computed  and  used 
in  the  empirical  scheme.  As  discussed  earlier,  the  MLD  pro- 
vides a  base  for  the  variables  DT5,  DT15,  DT50,  DT100,  DTf 
DZ,  DT3UZ  and  DZEUZ  to  define  their  values  in  the  predicted 
temperature  profile.  The  ether  variables  30T,  TBOT,  BOZ  and 
TBOZ,  besides  contributing  to  the  construction  of  the  pro- 
file, provide  a  check  to  maintain  continuity  with  the  MLD- 
determined  variables.  The  technique  used  in  this  thesis  to 
compute  MLD  caused  no  problems  in  the  determination  of  the 
thermal  profile;  however  a  comparison  test  of  computing 
mixed  layer  depth  by  the  five  different  methods,  defined 
earlier  in  Chapter  IV,  was  conducted  to  ascertain  which 
method   is    the    nicst   accurate. 

All  methods  viere  provided  wi-h  the  same  SST  values  from 
the  original  data  of  -he  ncrth-scuth  track,  and  compared  with 
the  actual  mix  =  d  layer  depth.  The  computed  depths  were  then 
correlated   to      determine    the    best   technique.  Results    from 

the  test  indicated  all  methods  were  surprisingly  equal.  The 
methed  used  was  chosen  because  of  the  higher  correlations 
between  variables.  However,  the  simplist  method,  which  was 
the  SST-derived  lixed  layer  depth,  with  a  correlation  coef- 
ficient   of    0.535,   could    have    been    used   with    equal    success. 
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v-       CONCLUSIONS    AND     RECOMMENDATIONS 

A.       CONCLUSIONS 

The  s€3  surface  temperature  reflects  many  of  the  physi- 
cal processes  that  govern  the  oceans.  Use  of  remote  sensing 
-co  map  the  surface  thermal  structure  of  the  world's  oceans 
accurately  and  rapidly  presently  exists,  but  the  determina- 
tion of  precise  subsurface  vertical  temperature  structure 
from  these  surface  signatures  is  unrealized.  Possibly  cnly 
in  selected  ocean  regions  can  insight  be  gained  into  deter- 
mining quantitative  relationships  for  directly  inferring  the 
subsurface  features  from  the  sea  surface  temperature.  In 
relatively  quiet  locations  with  uniformly  changing  proper- 
ties, such  as  in  the  Transitional  Zone  of  the  Northeast 
Pacific,  many  strong  correlations  can  be  observed  within  the 
thermal    structure. 

The  results  of  the  statistical  analysis  indicated  that 
the  strongest  relationships  within  the  vertical  temperature 
prcfile  cccurred  with  the  thermally-dependent  variables 
(TBOT,  DI,  TEUZ,  DTBUZ)  as  they  varied  with  SST.  The 
depth-related    variablis    (BLD,    30T,    DZ,    3UZ,    DZBUZ)     displayed 
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weaker  correlations  with  th€  surface  temperature,  but  trends 
still  could  be  observed  through  the  thermal  structure.  The 
use  of  scatter  diagrams  tc  assist  in  the  correlation  analy- 
sis cf  variables  was  necessary  to  establish  regional  bounda- 
ries and  the  best  fit  points  for  the  derivation  of 
regression   equations. 

An  important  region  within  the  vertical  temperature  pro- 
file is  the  thermocline.  In  practical  terms,  its  depth, 
strength  and  gradient  strcrgly  influence  the  propagation  cf 
underwater  sound.  Tne  prediction  of  this  structure  by  use 
of  only  sea  surface  temperature  would  greatly  enhance  ASW 
capability.  Normally,  within  the  Northeast  Pacific,  the 
correlation  between  SST  and  thermocline  depth  is  net  very 
high;  however,  in  the  Transition  Zone  with  uniformly  chang- 
ing structure  a  reasonable  predictability  can  be  obtained. 
The  thermocline  variables  (DT5,  DT15,  DT50,  DT103)  ail  pro- 
vided trends  frci  which  SST-related  values  were  computed  tc 
define  the  thermal  structure.  Overall,  the  mest  precisely 
predicted  region  of  the  temperature  profile  was  the  thermo- 
cline. The  depth  may  have  varied  due  to  internal  waves,  but 
the  negative  temperature  gradient  was  accurate  and  realistic 
when    compared   to    the    observed   structure. 
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B.        BECOMMENEATICNS 

In  future  research  cf  this  method,  an  expanded  sat  cf 
vertical  tsmpe rat ure  data  from  parallel  or  bisecting  -racks 
should  be  established  as  the  data  base  for  statistical  anal- 
ysis. The  addition  cf  an  east-west  track  of  data  would  pro- 
vide latitudinal  trends  which  could  also  be  used  in  the 
predicted  temperature  profiles.  As  in  any  data  aquisition, 
a  consider aticn  as  to  the  minimum  amount  cf  data  points  to 
determine  the  thermal  structure  statistically  should  be 
made.  Meanwhile,  a  number  of  no  less  than  20  profiles  ever 
a  distance  of  1000  km  with  spacing  between  points  not  to 
exceed   50    km    is    recommended. 

Other  suggestions  for  model  improvement  include  adding 
dynamical  and  thermodynamical  variables  such  as  wind  speed, 
direction  and  duration,  air  temperature,  radiation  flux,  and 
other  atmospheric  effects  *hich  could  alter  the  ocean  ther- 
mal structure.  Most  of  these  factors  can  be  determined  by 
remote  sensing;  data  from  other  satellite  sensors,  such  as 
an  altimeter,  SAS  and  Coastal  Zone  Color  Scanner  which  can 
measure   sea    surface    signatures,    should  also   be    included. 

lo  determine  the  limitation  of  accuracy  for  the  statis- 
tically-determined     profile      method,      more      stringent      tests 
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shculd  be  conducted  which  increase  the  spatial  and  time  fac- 
tors. In  certain  ocean  regions  one  could  possibly  maintain 
reasonable  accuracies  within  hundreds  of  kilometers  and  sev- 
eral weeks  duration  from  the  acquisition  of  thermal  data. 
Different  ccean  locations  should  be  analyzed  and  tested  to 
determine  the  pcssible  regions  where  the  thermal  profiles 
vary  in  an  orderly  fashion  which  can  be  easily  modeled. 
Areas  where  mescscale  features,  such  as  eddies  and  fronts, 
are  constantly  altering  the  water  column  by  advactive  pro- 
cesses  would    be    poor    candidates    for    this    method. 

Another  recommendation  is  to  test  in  other  seasons  of 
the  year  to  determine  if  the  statistical  approach  would 
apply.  The  difference  in  the  vertical  thermal  structure 
from  the  heating  to  cooling  season  is  so  great  that  possibly 
there  could  be  only  a  certain  window  of  time  that  allows  the 
sea    surface   temperature    to   define   the    subsurface    structure. 

Finally,  a  method  for  determining  the  occurrence  of 
large  temperature  inversions  in  the  Subarctic  and  Transition 
Regions  is  needed.  This  important  thermal  feature,  because 
of  its  variation  in  magnitude,  depth  and  location,  affects 
greatly  the  propagation  of  sound  by  forming  a  sound  channel 
immediately    below  the   thermccline.      If   this    structure    can    be 
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predicted  th<=n  the  full  potential  of  the  statistical  me-hca 
of  determining  subsurface  Thermal  structure  from  sea  surface 
temperature    wculd  be    realized. 
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2FPENDIX    A 

THERMAL    VARIABLES    VEESUS    LATITUDE/DISTANCE    ALONG 
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Figure      A.1.      SST  Versus    Latitude/Distance. 
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Figure      A.  2.      MLE  Versus    Latitude/Distance, 


1  19 


DISTANCE  (KM) 
100  300  300  400  S00  400  700  100  900  1000  1100 

■  II  I  II  Mi I.   .1.1.  I 


44  4)  43  41 

NORTH  LATITUDE  (DEO) 


Figure      A. 3.      D15   Versus    Latitude/Distance, 
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Figure      A.4.      CT15   Versus    Latitude/Distance, 
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Figure      A.  5.      ET50   Versus   Latitude/Distance. 
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Figure      A. 6.      ET100    Versus   Latit ade/Distanca 
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Figure      A.  11.      BUZ   Versus    Latitude/ Distance. 
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APPENDIX    B 

BIOMEDICAL   COMPUTES    PROGRAMS     (DIXON    AND    BROWN, 

1979) 

The  EICMED  statistical  computer  programs  were  utilized 
because  of  the  simple  and  flexible  programming  it  provided 
in  analyzing  the  thermal  data.  In  constructing  the  vertical 
temperature  profiles,  the  statistical  analysis  required  cor- 
relation and  regression  equation  computation  and  scatter 
diagram    plotting   of    the    thermal    variabiss. 

Ths    program    initially    selected      was    'Stepwise   Regression 

P2R.1       This      program    ccirputes      simple    statistics,         i.=., 

mean,    standard    deviation,    etc.,    correlation   coefficients    and 

linear      regression   equations      for    any      number   of      variables. 

The    statistics    used    in   the   computations   are    defined   as: 


•      mean 


.      standard    deviation  S   =        ^  (  Xj  -  X  )*/(  N  -  1  j 


—  vi 


-,& 


•      the    correlation    coefficient    between    two    variables 


r  - 
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In  simple  linear  regression  analysis,  where  there  is 
only  cna  independent  variable,  the  coefficients  are  deter- 
mined by  a  least  squares  method.  The  form  cf  ths  equations 
are: 

y  =  A  +  Bx   (where  A  and  B  are  the  linear  coefficients) 

The  EICMEE  computer  program  used  to  print  the  regional 
scatter  diagrams  is  titled  *Bivariate  (Scatter)  Plots 
P6D. '  Recressior  equations  for  regressing  variable  x  on  y 
and  y  on  x  can  be  requested  along  with  the  intersections  of 
the  regression  lines  with  the  axis  of  the  plot.  This  was 
useful  in  analyzing  the  data  zo  determine  the  best  fit  for 
the  various  reaiens. 
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JEPENDIX    C 

REGIONAL    SCATTER    DIAGRAMS    -    SST    VERSUS    THERMAL 

VARIABLES 
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Figure   C.8.   Begional  Scatter  Diagram  -  SST  Versus  DT 


141 


LEGEND 

a 

•  -  SU3ASCTIC 

a 

3 

S 

a  -  TRANSITIONAL 

a 

x- FRONT 

a.  SUBTROPICAL 

a 

OS 

a 

1- 

a 

9 

* 

• 

a 

a 

*s 

o* 

e 

-• 

13 

• 

• 

A 

Ss- 

S* 

• 

Ui 

A 

a 

UL 

• 

A 

Q 

A 

jM- 

A 

A.  * 

•. 

.AA 

U* 

A 

A 

a 

A 

i 

X 

M 

Oo. 

X 

M 

A 

»— — — — '  »  '     ■  ■■■   i    iiiiiw ii  '  ' 

14  !ft  18  20  22  24 

SEA  SURFACE  TEMPERATURE  (Q 


Figure      C.9.       Eegional  Scat-er    Diagram    -    SST    Versus   DZ 


142 


• 

X 

o 

©- 
T 

A 

i 

• 

•                                              A 

V 

•  • 

Ul 

A 

Z 

3? 

A 

X 

o  ' 

A 

• 

a         a 

IU 

A 

P 

A 

A 

Zo 

• 

<T 

S 

A 

u. 

o 

A 

2° 

A 

% 

o 

a 

LEGEND 

o 

•  - 

SUBARCTIC 

x 

A« 

TRANSITIONAL 

a 

£o 

X« 

FRONT 

*-  • 

Om 

SUBTROPICAL 

a 

o  ' 

a 

S. 

• 

a 

«  V4  14  II  80  n  2* 

SEA  SURFACE  TEMPERATURE  CO 


Figure   C.10.   Regional  Scatter  Diagram  -  S5T  Versus  BUZ 


143 


IU 

Z 

-j 

O 


UJ 

is 

u. 
O 

o 

t 

o 


< 

a. 

i 


•-* 


a   a 
a 


■xd  a 


AA 


A 


A  LEGEND 

•  -  SU3ABCTIC 


•  a  -  TRANSITIONAL 

•  x  »  FRONT 

#  a  -  SUBTROPICAL 


B  M  14  ft  20  32  M 

SEA  SURFACE  TEMPERATURE  (O 


Figurs   C.11.   Regional  Scatter  Diagram  -  SST  Versus  TBUZ 


144 


M 

3_ 


9 

O 

2 

O 

z 
< 

x 


3 


A 

A 
A 


A 


•  •  A 


a 

aa 
a 

a 


4V 

*  AD 


LEGEND 

•  -  su3Ancnc 

a  -  T2ANS1TI0NAL 

x  -  FRONT 

a  -  SUBTROPICAL 


■"->*■—■ — — —W— — ——H  |  I     —— — a—  ■■!■  IIMI  H  II  J      — — —         ■   II 

IS  14  16  l«  20  22  24 

SEA  SURFACE  TEMPERATURE  CO 


p- 


igure      C.12.      Begional    Scatter    Diagram   -    SST    Versus    DTBJZ 


145 


^  « 

LEGEND 

a 

2  ° 

•  -  SUBARCTIC 

a 

a  -  T2ANSIT10NAL 

a 

3 

x  -  FRONT 

SB 

a- SUBTROPICAL 

00 

►-S- 

9*" 

a 

z 

* 

A 

a 

OS 

£2 

& 

Ul 

a 

z 

• 

A 

a 

Ul 

ate 

A 

jg=- 

& 

u. 
a 

• 

A 

A        X 

X 

a 

k- 

°»  o 

1 

\     *  4 

N 

•    * 

3 

ca 

x 

h*o. 

aK 

• 

13 


1*  tS  20 

SEA  SURFACE  TEMPERATURE  CO 


32 


34 


Figure      C.  13.       Fsgicnal    Scatter    Diagram  -    55T    Versus    DZBUZ 


146 


APPENDIX    D 

REGIONAI    SCATTER    DIAGRAMS    -    EOT    VERSUS    THERMAL 

VARIABLES 
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REGIONAL    SCATTER    DIAGRAMS    -    BUZ    VERSUS    THERMA! 
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Figure   E.  1.   Regional  Scatter  Diagram  -  BUZ  Versus  DT5. 
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